<<<<<<< HEAD ======= >>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4 ======= >>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
<<<<<<< HEAD

Data

library(survival)# the package for durvival functions.
library(ggfortify) # Usefull for intepreting Survival functions plots.
## Loading required package: ggplot2
library(KMsurv)
Cdata<-read.table("C_data.txt",header = FALSE)
head(Cdata)
=======

Data

library(survival)# the package for durvival functions.
## Warning: package 'survival' was built under R version 4.0.3
library(ggfortify) # Usefull for intepreting Survival functions plots.
## Warning: package 'ggfortify' was built under R version 4.0.3
## Loading required package: ggplot2
Cdata<-read.table("C_data.txt",header = FALSE) #loading our data
head(Cdata) 
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##    V1 V2 V3 V4 V5 V6 V7 V8
## 1 204  1  4  0  1  2 66 50
## 2 374  1  5  0  0  2 66 43
## 3 184  1 13  1  1  2 57 43
## 4  60  1 16  1  0  2 67 40
## 5 449  1 21  1  1  2 62 25
## 6 224  0 22  0  0  0 75 94
<<<<<<< HEAD
colnames(Cdata) <-c("Patient_number", "Death_Censoring", "Time", "Treatment", "Sex", "Ascites", "Age", "Prothrombin")
=======
colnames(Cdata) <-c("Patient_number", "Death_Censoring", "Time", "Treatment", "Sex", "Ascites", "Age", "Prothrombin") #adding names
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
head(Cdata)
##   Patient_number Death_Censoring Time Treatment Sex Ascites Age Prothrombin
## 1            204               1    4         0   1       2  66          50
## 2            374               1    5         0   0       2  66          43
## 3            184               1   13         1   1       2  57          43
## 4             60               1   16         1   0       2  67          40
## 5            449               1   21         1   1       2  62          25
## 6            224               0   22         0   0       0  75          94
<<<<<<< HEAD

Opgave 1

Surv(Cdata$Time, Cdata$Death_Censoring)[1:10]
##  [1]  4   5  13  16  21  22+ 23  24+ 24+ 25+
sur.fit1 <- survfit(Surv(Time, Death_Censoring) ~ Treatment, data=Cdata)
=======

Opgave 1

sur.fit1 <- survfit(Surv(Time, Death_Censoring) ~ Treatment, data=Cdata)# constructing Survival function
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
names(sur.fit1)
##  [1] "n"         "time"      "n.risk"    "n.event"   "n.censor"  "surv"     
##  [7] "std.err"   "cumhaz"    "std.chaz"  "strata"    "type"      "logse"    
## [13] "conf.int"  "conf.type" "lower"     "upper"     "call"
<<<<<<< HEAD
plot(sur.fit1,xlab="Days",ylab="Overall survival probability")

autoplot(sur.fit1,xlab="Days",ylab="Overall survival probability")

=======
#plot(sur.fit1,xlab="Days",ylab="Overall survival probability") #a normal plot
autoplot(sur.fit1,xlab="Days",ylab="Overall survival probability") # a clear plot

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
# Log-rank test
Log_rank_test<-survdiff(Surv(Time, Death_Censoring) ~Treatment , data = Cdata)
Log_rank_test # results
## Call:
## survdiff(formula = Surv(Time, Death_Censoring) ~ Treatment, data = Cdata)
## 
##               N Observed Expected (O-E)^2/E (O-E)^2/V
## Treatment=0 251      142      149     0.355     0.728
## Treatment=1 237      150      143     0.371     0.728
## 
##  Chisq= 0.7  on 1 degrees of freedom, p= 0.4
<<<<<<< HEAD
1 - pchisq(Log_rank_test$chisq, length(Log_rank_test$n) - 1)#p-value
## [1] 0.3936019

Opgave 2

=======
1 - pchisq(Log_rank_test$chisq, length(Log_rank_test$n) - 1)#True p-value
## [1] 0.3936019

The log-rank test show no significans in Treatment method

Opgave 2

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
#Looking at Treatment
Cox.fit1 <- coxph(Sur1~ strata(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
summary(Cox.fit1)
## Call:
## coxph(formula = Sur1 ~ strata(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                       coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Sex)1      0.517727  1.678209  0.127220  4.070 4.71e-05 ***
## factor(Ascites)1  0.439503  1.551935  0.180240  2.438   0.0148 *  
## factor(Ascites)2  0.966760  2.629412  0.182675  5.292 1.21e-07 ***
## Age               0.048929  1.050145  0.006810  7.185 6.74e-13 ***
## Prothrombin      -0.012629  0.987450  0.002916 -4.331 1.49e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                  exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1        1.6782     0.5959    1.3078    2.1535
## factor(Ascites)1    1.5519     0.6444    1.0901    2.2095
## factor(Ascites)2    2.6294     0.3803    1.8381    3.7614
## Age                 1.0501     0.9522    1.0362    1.0643
## Prothrombin         0.9875     1.0127    0.9818    0.9931
## 
## Concordance= 0.692  (se = 0.018 )
## Likelihood ratio test= 126.1  on 5 df,   p=<2e-16
## Wald test            = 128.9  on 5 df,   p=<2e-16
## Score (logrank) test = 138.4  on 5 df,   p=<2e-16
source("andersen.R")
test1<-survfit(Cox.fit1)
test1
## Call: survfit(formula = Cox.fit1)
## 
##               n events median 0.95LCL 0.95UCL
## Treatment=0 251    142   1866    1374    2376
## Treatment=1 237    150   1588    1185    1978
andersen.plot(test1) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD

=======

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##   [1] 0.002803319 0.005646071 0.005646071 0.005646071 0.005646071 0.005646071
##   [7] 0.008523343 0.008523343 0.008523343 0.011437862 0.017335724 0.017335724
##  [13] 0.020345066 0.023395335 0.023395335 0.023395335 0.026499061 0.029643468
##  [19] 0.029643468 0.029643468 0.036035501 0.036035501 0.039307517 0.039307517
##  [25] 0.039307517 0.039307517 0.042649418 0.042649418 0.046063196 0.046063196
##  [31] 0.046063196 0.046063196 0.049541013 0.049541013 0.049541013 0.053049853
##  [37] 0.056566027 0.056566027 0.060164114 0.063845205 0.063845205 0.063845205
##  [43] 0.067587279 0.071356145 0.071356145 0.071356145 0.071356145 0.075181874
##  [49] 0.079202607 0.083279132 0.087385608 0.087385608 0.091568985 0.095774186
##  [55] 0.095774186 0.095774186 0.095774186 0.099991040 0.099991040 0.104263418
##  [61] 0.112937833 0.117394860 0.117394860 0.121898044 0.121898044 0.121898044
##  [67] 0.121898044 0.121898044 0.121898044 0.126483786 0.126483786 0.126483786
##  [73] 0.126483786 0.131091385 0.131091385 0.131091385 0.140583177 0.140583177
##  [79] 0.145384081 0.150193795 0.150193795 0.150193795 0.155024697 0.155024697
##  [85] 0.155024697 0.164843078 0.169809657 0.174879989 0.174879989 0.174879989
##  [91] 0.179983715 0.185148874 0.185148874 0.185148874 0.190343605 0.190343605
##  [97] 0.195754728 0.201215599 0.201215599 0.201215599 0.206791483 0.212387598
## [103] 0.212387598 0.212387598 0.212387598 0.218077148 0.218077148 0.223830989
## [109] 0.223830989 0.223830989 0.229677281 0.235643235 0.235643235 0.241738991
## [115] 0.241738991 0.241738991 0.247896302 0.247896302 0.254199298 0.254199298
## [121] 0.254199298 0.254199298 0.260596686 0.260596686 0.267079498 0.267079498
## [127] 0.273660243 0.273660243 0.273660243 0.280259522 0.280259522 0.280259522
## [133] 0.286921714 0.286921714 0.286921714 0.293635839 0.293635839 0.300422372
## [139] 0.300422372 0.300422372 0.300422372 0.300422372 0.300422372 0.307299696
## [145] 0.314189319 0.314189319 0.314189319 0.314189319 0.314189319 0.314189319
## [151] 0.314189319 0.321146855 0.328152411 0.335176681 0.335176681 0.342270564
## [157] 0.342270564 0.342270564 0.342270564 0.342270564 0.349485048 0.349485048
## [163] 0.356733711 0.356733711 0.356733711 0.364048567 0.364048567 0.364048567
## [169] 0.364048567 0.364048567 0.371491815 0.371491815 0.371491815 0.379164770
## [175] 0.379164770 0.386915944 0.386915944 0.394805438 0.402733066 0.402733066
## [181] 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066
## [187] 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066
## [193] 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066
## [199] 0.402733066 0.411236396 0.411236396 0.411236396 0.411236396 0.411236396
## [205] 0.411236396 0.411236396 0.411236396 0.419795452 0.437088202 0.437088202
## [211] 0.437088202 0.437088202 0.437088202 0.437088202 0.437088202 0.437088202
## [217] 0.437088202 0.437088202 0.437088202 0.437088202 0.437088202 0.446533320
## [223] 0.446533320 0.456102394 0.465779829 0.465779829 0.465779829 0.465779829
## [229] 0.475600891 0.485506021 0.495492038 0.495492038 0.495492038 0.495492038
## [235] 0.495492038 0.495492038 0.505788223 0.505788223 0.516238757 0.526874067
## [241] 0.537591327 0.537591327 0.537591327 0.548523094 0.548523094 0.559879460
## [247] 0.571271313 0.571271313 0.571271313 0.582703773 0.582703773 0.594325938
## [253] 0.594325938 0.594325938 0.594325938 0.594325938 0.594325938 0.594325938
## [259] 0.594325938 0.606511480 0.606511480 0.618765886 0.631204366 0.631204366
## [265] 0.631204366 0.631204366 0.643858883 0.643858883 0.643858883 0.643858883
## [271] 0.656704522 0.669715610 0.669715610 0.682993370 0.696337571 0.696337571
## [277] 0.696337571 0.696337571 0.696337571 0.696337571 0.710249295 0.710249295
## [283] 0.710249295 0.710249295 0.710249295 0.710249295 0.710249295 0.710249295
## [289] 0.724349453 0.724349453 0.724349453 0.724349453 0.724349453 0.738577046
## [295] 0.738577046 0.753446696 0.753446696 0.753446696 0.753446696 0.753446696
## [301] 0.753446696 0.753446696 0.753446696 0.768898259 0.768898259 0.768898259
## [307] 0.768898259 0.768898259 0.768898259 0.768898259 0.768898259 0.785192921
## [313] 0.785192921 0.785192921 0.785192921 0.801765910 0.818684550 0.818684550
## [319] 0.818684550 0.836023320 0.836023320 0.853639281 0.871369433 0.889301550
## [325] 0.889301550 0.889301550 0.907645633 0.907645633 0.926439496 0.926439496
## [331] 0.926439496 0.926439496 0.926439496 0.947256706 0.947256706 0.947256706
## [337] 0.947256706 0.947256706 0.947256706 0.947256706 0.947256706 0.947256706
## [343] 0.947256706 0.947256706 0.969892125 0.969892125 0.969892125 0.969892125
## [349] 0.969892125 0.969892125 0.969892125 0.969892125 0.969892125 0.994528565
## [355] 0.994528565 0.994528565 0.994528565 0.994528565 0.994528565 0.994528565
## [361] 1.020675491 1.020675491 1.020675491 1.020675491 1.020675491 1.049912169
## [367] 1.079375654 1.079375654 1.079375654 1.079375654 1.079375654 1.079375654
## [373] 1.079375654 1.079375654 1.079375654 1.146811358 1.146811358 1.181835055
## [379] 1.218516764 1.218516764 1.256951458 1.298382879 1.298382879 1.341388730
## [385] 1.341388730 1.341388730 1.386772376 1.433078818 1.480546087 1.480546087
## [391] 1.480546087 1.529820305 1.529820305 1.529820305 1.529820305 1.529820305
## [397] 1.529820305 1.529820305 1.529820305 1.529820305 1.529820305 1.529820305
## [403] 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801
## [409] 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801
## [415] 1.588748801 1.588748801 1.588748801 1.706570781 1.851428550 1.851428550
## [421] 1.851428550 1.851428550 2.082565012 2.082565012 2.082565012 2.082565012
## [427] 2.082565012 2.082565012 2.082565012 2.082565012 2.082565012 2.082565012
## [433] 0.000000000 0.000000000 0.003250690 0.006538576 0.009864885 0.009864885
## [439] 0.009864885 0.009864885 0.009864885 0.013277201 0.013277201 0.013277201
## [445] 0.013277201 0.013277201 0.020207038 0.020207038 0.020207038 0.020207038
## [451] 0.020207038 0.023702448 0.027228455 0.030806238 0.034435535 0.038074158
## [457] 0.038074158 0.038074158 0.038074158 0.041732456 0.041732456 0.045403818
## [463] 0.049098253 0.049098253 0.049098253 0.052808653 0.056530257 0.056530257
## [469] 0.056530257 0.056530257 0.056530257 0.056530257 0.060256972 0.060256972
## [475] 0.060256972 0.060256972 0.060256972 0.060256972 0.060256972 0.071604831
## [481] 0.071604831 0.071604831 0.071604831 0.071604831 0.071604831 0.071604831
## [487] 0.079371547 0.083414992 0.087462234 0.087462234 0.099766915 0.099766915
## [493] 0.099766915 0.099766915 0.099766915 0.099766915 0.103915026 0.108113806
## [499] 0.108113806 0.108113806 0.108113806 0.108113806 0.112339438 0.112339438
## [505] 0.116587624 0.116587624 0.120853383 0.125300922 0.134296069 0.138899626
## [511] 0.138899626 0.143534473 0.143534473 0.148228377 0.148228377 0.152996332
## [517] 0.152996332 0.152996332 0.152996332 0.152996332 0.152996332 0.157794374
## [523] 0.162621192 0.162621192 0.162621192 0.167501229 0.172394281 0.177313034
## [529] 0.177313034 0.177313034 0.182249823 0.187236669 0.187236669 0.187236669
## [535] 0.192250769 0.197293558 0.202362380 0.207474451 0.207474451 0.207474451
## [541] 0.207474451 0.207474451 0.207474451 0.207474451 0.212634116 0.212634116
## [547] 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116
## [553] 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116
## [559] 0.212634116 0.218014502 0.218014502 0.218014502 0.223448922 0.223448922
## [565] 0.223448922 0.223448922 0.234510155 0.234510155 0.234510155 0.234510155
## [571] 0.240106474 0.245734315 0.251402356 0.251402356 0.251402356 0.251402356
## [577] 0.251402356 0.251402356 0.262881729 0.268674343 0.274480410 0.280323394
## [583] 0.286199096 0.286199096 0.286199096 0.286199096 0.286199096 0.286199096
## [589] 0.292129315 0.298073929 0.304107082 0.304107082 0.304107082 0.310162759
## [595] 0.310162759 0.316318573 0.322505090 0.322505090 0.322505090 0.322505090
## [601] 0.322505090 0.322505090 0.322505090 0.329083016 0.329083016 0.329083016
## [607] 0.335775894 0.335775894 0.342588674 0.342588674 0.342588674 0.342588674
## [613] 0.349410916 0.349410916 0.356306234 0.356306234 0.363250084 0.370238570
## [619] 0.377324446 0.377324446 0.384440778 0.391586591 0.391586591 0.398786655
## [625] 0.406261565 0.413790884 0.421370627 0.428975324 0.428975324 0.436859628
## [631] 0.444804484 0.444804484 0.452830595 0.460901444 0.469010410 0.477281466
## [637] 0.485627288 0.494045343 0.494045343 0.494045343 0.494045343 0.502797105
## [643] 0.511776939 0.511776939 0.511776939 0.511776939 0.520903087 0.530109827
## [649] 0.530109827 0.539428699 0.548789215 0.558218793 0.558218793 0.558218793
## [655] 0.558218793 0.567719358 0.567719358 0.567719358 0.577406531 0.587144719
## [661] 0.587144719 0.587144719 0.587144719 0.587144719 0.587144719 0.587144719
## [667] 0.587144719 0.597190566 0.597190566 0.607392500 0.607392500 0.607392500
## [673] 0.607392500 0.607392500 0.618050787 0.618050787 0.618050787 0.618050787
## [679] 0.618050787 0.618050787 0.628915165 0.628915165 0.628915165 0.628915165
## [685] 0.628915165 0.628915165 0.640015400 0.651237007 0.662783323 0.674765510
## [691] 0.686841825 0.686841825 0.699071671 0.699071671 0.699071671 0.699071671
## [697] 0.711603645 0.724294283 0.724294283 0.737229093 0.750560072 0.764056687
## [703] 0.764056687 0.764056687 0.777823629 0.777823629 0.777823629 0.777823629
## [709] 0.777823629 0.777823629 0.792009479 0.806448636 0.806448636 0.806448636
## [715] 0.821296645 0.821296645 0.836526845 0.852054132 0.867775447 0.883773968
## [721] 0.883773968 0.899977247 0.916477011 0.933106732 0.950032324 0.950032324
## [727] 0.950032324 0.950032324 0.950032324 0.950032324 0.950032324 0.950032324
## [733] 0.950032324 0.950032324 0.950032324 0.950032324 0.950032324 0.968763551
## [739] 0.987794194 1.007279790 1.007279790 1.007279790 1.027323522 1.027323522
## [745] 1.027323522 1.027323522 1.048040964 1.048040964 1.048040964 1.048040964
## [751] 1.048040964 1.048040964 1.069921181 1.069921181 1.069921181 1.069921181
## [757] 1.069921181 1.069921181 1.069921181 1.069921181 1.069921181 1.094575178
## [763] 1.119831131 1.119831131 1.148143865 1.148143865 1.148143865 1.177604111
## [769] 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111
## [775] 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111
## [781] 1.209863317 1.243221810 1.243221810 1.243221810 1.243221810 1.243221810
## [787] 1.277550368 1.277550368 1.313206852 1.313206852 1.313206852 1.351565798
## [793] 1.351565798 1.351565798 1.351565798 1.351565798 1.351565798 1.351565798
## [799] 1.351565798 1.351565798 1.351565798 1.351565798 1.351565798 1.392826090
## [805] 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090
## [811] 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090
## [817] 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090
## [823] 1.446508714 1.446508714 1.446508714 1.446508714 1.446508714 1.446508714
## [829] 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976
## [835] 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976
## [841] 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976
## [847] 1.669782450 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895
## [853] 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895
## [859] 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895
res1<-residuals(Cox.fit1,type="deviance")
hist(res1)
<<<<<<< HEAD

boxplot(res1~Cdata$Death_Censoring)

=======

boxplot(res1~Cdata$Death_Censoring)

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
#Looking at sex
Cox.fit2 <- coxph(Sur1~ factor(Treatment)+ strata(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
summary(Cox.fit2)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + strata(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                         coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Treatment)1  0.054861  1.056393  0.118279  0.464   0.6428    
## factor(Ascites)1    0.455062  1.576270  0.179616  2.534   0.0113 *  
## factor(Ascites)2    0.977054  2.656619  0.182538  5.353 8.67e-08 ***
## Age                 0.048919  1.050135  0.006785  7.209 5.62e-13 ***
## Prothrombin        -0.012850  0.987233  0.002937 -4.376 1.21e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                    exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1    1.0564     0.9466    0.8378    1.3320
## factor(Ascites)1      1.5763     0.6344    1.1085    2.2414
## factor(Ascites)2      2.6566     0.3764    1.8576    3.7993
## Age                   1.0501     0.9523    1.0363    1.0642
## Prothrombin           0.9872     1.0129    0.9816    0.9929
## 
## Concordance= 0.703  (se = 0.018 )
## Likelihood ratio test= 125  on 5 df,   p=<2e-16
## Wald test            = 129.3  on 5 df,   p=<2e-16
## Score (logrank) test = 140.6  on 5 df,   p=<2e-16
source("andersen.R")
test2<-survfit(Cox.fit2)
test2
## Call: survfit(formula = Cox.fit2)
## 
##         n events median 0.95LCL 0.95UCL
## Sex=0 198    111   2455    1950    2994
## Sex=1 290    181   1373    1141    1729
andersen.plot(test2) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD

=======

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##   [1] 0.000000000 0.003345221 0.003345221 0.006739532 0.006739532 0.006739532
##   [7] 0.010207978 0.010207978 0.010207978 0.010207978 0.010207978 0.010207978
##  [13] 0.013778983 0.013778983 0.017427731 0.017427731 0.017427731 0.017427731
##  [19] 0.017427731 0.021084701 0.024789265 0.028503991 0.032298009 0.032298009
##  [25] 0.032298009 0.032298009 0.032298009 0.036167088 0.036167088 0.036167088
##  [31] 0.040056452 0.040056452 0.043970676 0.043970676 0.043970676 0.047938128
##  [37] 0.047938128 0.047938128 0.051917872 0.051917872 0.056005352 0.056005352
##  [43] 0.056005352 0.056005352 0.056005352 0.056005352 0.056005352 0.060109904
##  [49] 0.060109904 0.064231724 0.064231724 0.064231724 0.068394312 0.072585660
##  [55] 0.072585660 0.076812794 0.076812794 0.076812794 0.089582682 0.093887342
##  [61] 0.098227289 0.098227289 0.098227289 0.102626326 0.102626326 0.102626326
##  [67] 0.102626326 0.102626326 0.102626326 0.107108781 0.111612664 0.111612664
##  [73] 0.111612664 0.116159004 0.116159004 0.116159004 0.130114809 0.130114809
##  [79] 0.134930460 0.139757868 0.139757868 0.139757868 0.139757868 0.139757868
##  [85] 0.139757868 0.139757868 0.139757868 0.139757868 0.139757868 0.139757868
##  [91] 0.144647340 0.144647340 0.144647340 0.144647340 0.144647340 0.149551589
##  [97] 0.154481057 0.154481057 0.154481057 0.159466378 0.159466378 0.159466378
## [103] 0.164489861 0.169554440 0.174655924 0.179818794 0.179818794 0.185011809
## [109] 0.185011809 0.185011809 0.185011809 0.185011809 0.185011809 0.185011809
## [115] 0.185011809 0.185011809 0.185011809 0.185011809 0.190304939 0.190304939
## [121] 0.190304939 0.190304939 0.190304939 0.190304939 0.190304939 0.190304939
## [127] 0.190304939 0.195726388 0.195726388 0.195726388 0.201225033 0.201225033
## [133] 0.201225033 0.201225033 0.206830590 0.206830590 0.206830590 0.206830590
## [139] 0.206830590 0.206830590 0.212498753 0.212498753 0.212498753 0.218251218
## [145] 0.224014976 0.224014976 0.229841369 0.235707790 0.241593267 0.241593267
## [151] 0.247531873 0.253548931 0.253548931 0.253548931 0.253548931 0.253548931
## [157] 0.253548931 0.253548931 0.259625289 0.259625289 0.265733659 0.271874296
## [163] 0.271874296 0.271874296 0.278161128 0.278161128 0.278161128 0.278161128
## [169] 0.278161128 0.278161128 0.278161128 0.278161128 0.278161128 0.278161128
## [175] 0.278161128 0.278161128 0.284764700 0.291380525 0.291380525 0.291380525
## [181] 0.291380525 0.291380525 0.291380525 0.291380525 0.291380525 0.291380525
## [187] 0.291380525 0.291380525 0.298197782 0.305053024 0.305053024 0.312181122
## [193] 0.312181122 0.319696199 0.327281618 0.334901911 0.334901911 0.334901911
## [199] 0.334901911 0.342926589 0.350976673 0.359089262 0.367255347 0.367255347
## [205] 0.367255347 0.375658814 0.375658814 0.384589807 0.384589807 0.384589807
## [211] 0.384589807 0.384589807 0.384589807 0.384589807 0.394043327 0.394043327
## [217] 0.394043327 0.394043327 0.403874432 0.403874432 0.403874432 0.403874432
## [223] 0.403874432 0.413933604 0.413933604 0.413933604 0.413933604 0.413933604
## [229] 0.413933604 0.413933604 0.424213502 0.424213502 0.424213502 0.424213502
## [235] 0.424213502 0.424213502 0.424213502 0.424213502 0.434650119 0.445333746
## [241] 0.445333746 0.445333746 0.456279726 0.456279726 0.456279726 0.456279726
## [247] 0.467431862 0.467431862 0.467431862 0.467431862 0.467431862 0.478635809
## [253] 0.478635809 0.478635809 0.478635809 0.478635809 0.490022711 0.502016289
## [259] 0.502016289 0.514141966 0.514141966 0.514141966 0.514141966 0.514141966
## [265] 0.514141966 0.514141966 0.514141966 0.514141966 0.514141966 0.514141966
## [271] 0.514141966 0.514141966 0.514141966 0.526358492 0.538649435 0.538649435
## [277] 0.538649435 0.538649435 0.538649435 0.551125907 0.551125907 0.551125907
## [283] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [289] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [295] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [301] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [307] 0.577806544 0.577806544 0.577806544 0.577806544 0.577806544 0.592307003
## [313] 0.592307003 0.592307003 0.607040920 0.607040920 0.607040920 0.607040920
## [319] 0.607040920 0.607040920 0.607040920 0.622298172 0.637667908 0.653237280
## [325] 0.653237280 0.653237280 0.653237280 0.653237280 0.669499364 0.669499364
## [331] 0.687825832 0.687825832 0.708460386 0.729552770 0.729552770 0.750997675
## [337] 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675
## [343] 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675
## [349] 0.776640900 0.776640900 0.776640900 0.776640900 0.776640900 0.776640900
## [355] 0.803716523 0.803716523 0.803716523 0.803716523 0.803716523 0.833312019
## [361] 0.833312019 0.833312019 0.833312019 0.833312019 0.833312019 0.865278809
## [367] 0.865278809 0.865278809 0.865278809 0.865278809 0.865278809 0.899320721
## [373] 0.899320721 0.899320721 0.899320721 0.937040991 0.937040991 0.937040991
## [379] 0.937040991 0.937040991 0.975922999 0.975922999 0.975922999 1.018994173
## [385] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [391] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [397] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [403] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [409] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [415] 1.118012656 1.118012656 1.118012656 1.234982876 1.366600814 1.366600814
## [421] 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814
## [427] 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814
## [433] 0.002757535 0.002757535 0.005545448 0.005545448 0.008356406 0.008356406
## [439] 0.008356406 0.008356406 0.008356406 0.014074454 0.019861505 0.019861505
## [445] 0.019861505 0.022812378 0.025800620 0.025800620 0.028806331 0.031842293
## [451] 0.031842293 0.031842293 0.038015818 0.038015818 0.041181305 0.044352707
## [457] 0.044352707 0.044352707 0.047557766 0.047557766 0.050815624 0.054110208
## [463] 0.054110208 0.054110208 0.054110208 0.057430200 0.060757638 0.060757638
## [469] 0.064088457 0.064088457 0.064088457 0.067489583 0.067489583 0.067489583
## [475] 0.070931692 0.074391656 0.074391656 0.074391656 0.074391656 0.085107091
## [481] 0.088790926 0.088790926 0.092511874 0.092511874 0.092511874 0.092511874
## [487] 0.100137950 0.100137950 0.104080641 0.108053817 0.108053817 0.108053817
## [493] 0.112066311 0.116157297 0.116157297 0.116157297 0.120291324 0.124467166
## [499] 0.124467166 0.124467166 0.124467166 0.124467166 0.124467166 0.124467166
## [505] 0.128676372 0.128676372 0.132899856 0.137320430 0.141775934 0.146266732
## [511] 0.146266732 0.150782295 0.150782295 0.155326791 0.159929717 0.164567937
## [517] 0.164567937 0.173917460 0.178633284 0.183422929 0.183422929 0.188244562
## [523] 0.193090448 0.197980529 0.197980529 0.202928139 0.212933611 0.212933611
## [529] 0.212933611 0.218071708 0.223290566 0.223290566 0.228556180 0.233835838
## [535] 0.233835838 0.233835838 0.233835838 0.239180556 0.239180556 0.239180556
## [541] 0.239180556 0.239180556 0.244594070 0.250089607 0.255672777 0.261327717
## [547] 0.261327717 0.261327717 0.267073765 0.267073765 0.267073765 0.267073765
## [553] 0.267073765 0.267073765 0.272992839 0.272992839 0.278969209 0.278969209
## [559] 0.285010810 0.285010810 0.285010810 0.291064642 0.291064642 0.291064642
## [565] 0.297160780 0.297160780 0.303291099 0.309467267 0.309467267 0.315691443
## [571] 0.321969831 0.328280967 0.328280967 0.328280967 0.328280967 0.328280967
## [577] 0.328280967 0.328280967 0.334634071 0.334634071 0.334634071 0.341030821
## [583] 0.341030821 0.341030821 0.347459936 0.353901273 0.353901273 0.360380309
## [589] 0.366939999 0.373514305 0.373514305 0.373514305 0.373514305 0.373514305
## [595] 0.380178409 0.386886460 0.386886460 0.393624672 0.393624672 0.393624672
## [601] 0.393624672 0.393624672 0.400609212 0.407782608 0.407782608 0.415046290
## [607] 0.422365126 0.429802940 0.429802940 0.429802940 0.437342053 0.437342053
## [613] 0.444947615 0.444947615 0.452630371 0.452630371 0.460363210 0.468142839
## [619] 0.476022164 0.476022164 0.476022164 0.476022164 0.476022164 0.476022164
## [625] 0.483994395 0.483994395 0.483994395 0.483994395 0.483994395 0.492018105
## [631] 0.500093597 0.500093597 0.500093597 0.500093597 0.500093597 0.508238655
## [637] 0.516444093 0.516444093 0.516444093 0.516444093 0.533070772 0.541536117
## [643] 0.550180925 0.550180925 0.550180925 0.550180925 0.550180925 0.558938400
## [649] 0.558938400 0.567780493 0.567780493 0.576653625 0.576653625 0.585578905
## [655] 0.585578905 0.594577246 0.603650668 0.603650668 0.612858741 0.622104754
## [661] 0.631485851 0.640926781 0.640926781 0.640926781 0.640926781 0.640926781
## [667] 0.640926781 0.650667580 0.660529455 0.670502698 0.670502698 0.670502698
## [673] 0.680756087 0.680756087 0.680756087 0.691163264 0.691163264 0.701898395
## [679] 0.701898395 0.701898395 0.712689313 0.723583018 0.723583018 0.723583018
## [685] 0.723583018 0.723583018 0.734995825 0.746515007 0.746515007 0.746515007
## [691] 0.758324270 0.758324270 0.770255405 0.782427450 0.794741332 0.794741332
## [697] 0.807220965 0.819830004 0.832639334 0.845602503 0.858900126 0.872334656
## [703] 0.885990471 0.899793250 0.913836031 0.913836031 0.913836031 0.913836031
## [709] 0.913836031 0.913836031 0.928596975 0.928596975 0.943585879 0.943585879
## [715] 0.943585879 0.943585879 0.958987437 0.974640083 0.990455404 1.006504742
## [721] 1.022724644 1.039076112 1.055676826 1.072385926 1.089344418 1.106801640
## [727] 1.106801640 1.125021243 1.125021243 1.125021243 1.125021243 1.125021243
## [733] 1.125021243 1.125021243 1.125021243 1.144058979 1.144058979 1.163731681
## [739] 1.163731681 1.183675937 1.183675937 1.183675937 1.204101833 1.204101833
## [745] 1.204101833 1.204101833 1.204101833 1.225177632 1.246693520 1.246693520
## [751] 1.246693520 1.269109977 1.291888682 1.291888682 1.291888682 1.291888682
## [757] 1.291888682 1.291888682 1.316739239 1.316739239 1.316739239 1.342240244
## [763] 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244
## [769] 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244
## [775] 1.342240244 1.342240244 1.369335341 1.369335341 1.369335341 1.369335341
## [781] 1.369335341 1.397623227 1.397623227 1.397623227 1.397623227 1.426808943
## [787] 1.426808943 1.426808943 1.457454480 1.457454480 1.457454480 1.457454480
## [793] 1.488957202 1.488957202 1.488957202 1.488957202 1.488957202 1.488957202
## [799] 1.523228454 1.523228454 1.523228454 1.523228454 1.523228454 1.523228454
## [805] 1.523228454 1.523228454 1.523228454 1.560899224 1.560899224 1.600963914
## [811] 1.642737790 1.642737790 1.642737790 1.687062381 1.687062381 1.687062381
## [817] 1.687062381 1.687062381 1.735525345 1.784811988 1.835120582 1.835120582
## [823] 1.887003831 1.940567935 1.940567935 1.940567935 1.940567935 1.940567935
## [829] 2.006758591 2.006758591 2.006758591 2.006758591 2.006758591 2.006758591
## [835] 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732
## [841] 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732
## [847] 2.087886732 2.247980143 2.247980143 2.247980143 2.247980143 2.247980143
## [853] 2.247980143 2.247980143 2.599737118 2.599737118 2.599737118 2.599737118
## [859] 2.599737118 2.599737118 2.599737118 2.599737118 2.599737118 2.599737118
res2<-residuals(Cox.fit2,type="deviance")
hist(res2)
<<<<<<< HEAD

boxplot(res2~Cdata$Death_Censoring)

=======

boxplot(res2~Cdata$Death_Censoring)

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
# Looking at Ascites
Cox.fit3 <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ strata(Ascites)+ Age+ Prothrombin, data = Cdata)
summary(Cox.fit3)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + strata(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                         coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Treatment)1  0.092631  1.097057  0.118711  0.780  0.43521    
## factor(Sex)1        0.481967  1.619256  0.126547  3.809  0.00014 ***
## Age                 0.048120  1.049297  0.006795  7.082 1.42e-12 ***
## Prothrombin        -0.012599  0.987480  0.002929 -4.301 1.70e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                    exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1    1.0971     0.9115    0.8693    1.3844
## factor(Sex)1          1.6193     0.6176    1.2636    2.0751
## Age                   1.0493     0.9530    1.0354    1.0634
## Prothrombin           0.9875     1.0127    0.9818    0.9932
## 
## Concordance= 0.646  (se = 0.02 )
## Likelihood ratio test= 75.12  on 4 df,   p=2e-15
## Wald test            = 70.1  on 4 df,   p=2e-14
## Score (logrank) test = 70.43  on 4 df,   p=2e-14
source("andersen.R")
test3<-survfit(Cox.fit3)
test3
## Call: survfit(formula = Cox.fit3)
## 
##             n events median 0.95LCL 0.95UCL
## Ascites=0 386    211   1979    1618    2277
## Ascites=1  54     39   1194     802    2187
## Ascites=2  48     42    686     224    1529
andersen.plot(test3) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD

=======

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##    [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.002422910 0.002422910
##    [7] 0.002422910 0.002422910 0.007342705 0.007342705 0.009818343 0.012299473
##   [13] 0.012299473 0.017311969 0.019844167 0.022381833 0.022381833 0.022381833
##   [19] 0.024936882 0.024936882 0.027499712 0.030075870 0.030075870 0.032669946
##   [25] 0.035270579 0.037874232 0.037874232 0.037874232 0.040499620 0.043138885
##   [31] 0.043138885 0.045803991 0.045803991 0.045803991 0.045803991 0.045803991
##   [37] 0.051234459 0.053959659 0.053959659 0.056709190 0.056709190 0.056709190
##   [43] 0.059476974 0.062246981 0.067836269 0.067836269 0.067836269 0.070640906
##   [49] 0.070640906 0.073473549 0.076333898 0.076333898 0.076333898 0.076333898
##   [55] 0.079228569 0.082130814 0.082130814 0.085048413 0.087975882 0.087975882
##   [61] 0.090949553 0.096957268 0.099985028 0.103028855 0.109139889 0.109139889
##   [67] 0.109139889 0.112228053 0.112228053 0.118449962 0.121585833 0.121585833
##   [73] 0.124740883 0.127910783 0.131086222 0.131086222 0.134295773 0.137512089
##   [79] 0.140741665 0.140741665 0.140741665 0.143980878 0.147245615 0.150524681
##   [85] 0.153811349 0.153811349 0.157132578 0.160467711 0.167198221 0.167198221
##   [91] 0.170591826 0.170591826 0.170591826 0.170591826 0.174035043 0.174035043
##   [97] 0.174035043 0.174035043 0.174035043 0.174035043 0.174035043 0.177557058
##  [103] 0.177557058 0.181106991 0.181106991 0.184688635 0.188276315 0.188276315
##  [109] 0.188276315 0.191893457 0.191893457 0.195545936 0.195545936 0.199220479
##  [115] 0.202918422 0.202918422 0.206640238 0.210388559 0.214153663 0.217940133
##  [121] 0.217940133 0.217940133 0.221759485 0.225583152 0.225583152 0.233296347
##  [127] 0.237181790 0.241074706 0.244987967 0.248918660 0.252879175 0.256857068
##  [133] 0.260841555 0.260841555 0.260841555 0.264850422 0.268867250 0.272931469
##  [139] 0.272931469 0.277008479 0.281097733 0.285209776 0.289338169 0.293505342
##  [145] 0.293505342 0.293505342 0.293505342 0.293505342 0.297817238 0.297817238
##  [151] 0.302187933 0.306586088 0.311044871 0.315508708 0.319986281 0.319986281
##  [157] 0.319986281 0.319986281 0.324531278 0.324531278 0.324531278 0.329130385
##  [163] 0.333778961 0.333778961 0.338472789 0.343182554 0.347921280 0.352686126
##  [169] 0.357475797 0.362277850 0.367093534 0.371936042 0.376814323 0.381700710
##  [175] 0.386607441 0.386607441 0.391531514 0.396487204 0.396487204 0.401486385
##  [181] 0.406493245 0.406493245 0.411533108 0.411533108 0.411533108 0.411533108
##  [187] 0.416642840 0.421783672 0.421783672 0.427027876 0.432287963 0.437574903
##  [193] 0.437574903 0.442910860 0.442910860 0.453714946 0.459186349 0.459186349
##  [199] 0.464725859 0.464725859 0.470285325 0.475873942 0.481490529 0.481490529
##  [205] 0.481490529 0.481490529 0.481490529 0.487209959 0.492988895 0.498820289
##  [211] 0.498820289 0.504785833 0.504785833 0.510816417 0.510816417 0.516962041
##  [217] 0.523119100 0.523119100 0.529304753 0.535538855 0.535538855 0.541833971
##  [223] 0.541833971 0.541833971 0.548306163 0.548306163 0.554827461 0.561383167
##  [229] 0.567992686 0.567992686 0.574625244 0.581315779 0.581315779 0.588074133
##  [235] 0.594886993 0.601782986 0.601782986 0.608740305 0.615751400 0.622848256
##  [241] 0.622848256 0.629999870 0.637258708 0.644540098 0.644540098 0.644540098
##  [247] 0.644540098 0.652030508 0.652030508 0.659604587 0.659604587 0.667301465
##  [253] 0.667301465 0.675122748 0.683036373 0.691009874 0.699068259 0.707188520
##  [259] 0.715354708 0.723609767 0.731903802 0.740285562 0.748838785 0.748838785
##  [265] 0.757639804 0.757639804 0.757639804 0.757639804 0.757639804 0.757639804
##  [271] 0.757639804 0.757639804 0.757639804 0.757639804 0.766967351 0.776382646
##  [277] 0.785932909 0.785932909 0.785932909 0.795662298 0.805630556 0.805630556
##  [283] 0.805630556 0.815715614 0.825948817 0.836324102 0.836324102 0.836324102
##  [289] 0.847059561 0.857909699 0.868808764 0.879793185 0.879793185 0.879793185
##  [295] 0.891329087 0.891329087 0.903058135 0.903058135 0.914948511 0.926971810
##  [301] 0.926971810 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300
##  [307] 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300
##  [313] 0.939203300 0.939203300 0.952862217 0.966757789 0.966757789 0.966757789
##  [319] 0.966757789 0.981074755 0.995864031 0.995864031 1.010953971 1.010953971
##  [325] 1.010953971 1.026685796 1.042530764 1.042530764 1.042530764 1.042530764
##  [331] 1.042530764 1.059580899 1.076967045 1.076967045 1.076967045 1.076967045
##  [337] 1.076967045 1.076967045 1.076967045 1.076967045 1.116142086 1.116142086
##  [343] 1.136740519 1.136740519 1.136740519 1.158208569 1.180687663 1.180687663
##  [349] 1.204215922 1.204215922 1.204215922 1.229160124 1.254409050 1.280031424
##  [355] 1.280031424 1.306648637 1.333899572 1.333899572 1.333899572 1.333899572
##  [361] 1.333899572 1.333899572 1.333899572 1.333899572 1.333899572 1.333899572
##  [367] 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164
##  [373] 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164
##  [379] 1.430615609 1.498616602 1.498616602 1.571643693 1.655346574 1.655346574
##  [385] 1.655346574 1.655346574 1.780452655 1.780452655 1.780452655 1.780452655
##  [391] 1.780452655 1.780452655 1.780452655 1.780452655 1.780452655 0.000000000
##  [397] 0.012819623 0.012819623 0.012819623 0.026034998 0.039566407 0.039566407
##  [403] 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500
##  [409] 0.053402500 0.068122762 0.068122762 0.068122762 0.068122762 0.068122762
##  [415] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
##  [421] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
##  [427] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.100386010
##  [433] 0.100386010 0.100386010 0.100386010 0.117547315 0.134850559 0.134850559
##  [439] 0.134850559 0.134850559 0.153076521 0.171646666 0.171646666 0.171646666
##  [445] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
##  [451] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
##  [457] 0.191372504 0.191372504 0.191372504 0.191372504 0.191372504 0.211774184
##  [463] 0.211774184 0.211774184 0.211774184 0.211774184 0.211774184 0.211774184
##  [469] 0.232795610 0.232795610 0.232795610 0.232795610 0.232795610 0.232795610
##  [475] 0.254560972 0.276919329 0.276919329 0.276919329 0.276919329 0.276919329
##  [481] 0.300661572 0.300661572 0.300661572 0.300661572 0.300661572 0.300661572
##  [487] 0.300661572 0.300661572 0.324925481 0.324925481 0.324925481 0.324925481
##  [493] 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784
##  [499] 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784 0.378294064
##  [505] 0.378294064 0.378294064 0.378294064 0.378294064 0.406697290 0.406697290
##  [511] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
##  [517] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
##  [523] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
##  [529] 0.406697290 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806
##  [535] 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806
##  [541] 0.436020806 0.436020806 0.466837198 0.466837198 0.466837198 0.466837198
##  [547] 0.466837198 0.466837198 0.466837198 0.466837198 0.466837198 0.499883881
##  [553] 0.499883881 0.499883881 0.499883881 0.534298268 0.534298268 0.534298268
##  [559] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
##  [565] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
##  [571] 0.569580748 0.569580748 0.569580748 0.569580748 0.569580748 0.607417459
##  [577] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
##  [583] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
##  [589] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
##  [595] 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010
##  [601] 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010 0.756226472
##  [607] 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
##  [613] 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
##  [619] 0.756226472 0.756226472 0.818519693 0.818519693 0.818519693 0.818519693
##  [625] 0.889304890 0.889304890 0.889304890 0.889304890 0.889304890 0.889304890
##  [631] 0.889304890 0.968545149 0.968545149 0.968545149 0.968545149 1.060605655
##  [637] 1.060605655 1.060605655 1.060605655 1.060605655 1.060605655 1.060605655
##  [643] 1.060605655 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
##  [649] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
##  [655] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
##  [661] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
##  [667] 1.162895503 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [673] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [679] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [685] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [691] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [697] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [703] 1.288056766 1.288056766 1.288056766 1.472408461 1.472408461 1.472408461
##  [709] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [715] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [721] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [727] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [733] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [739] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877
##  [745] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877
##  [751] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.986477110
##  [757] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
##  [763] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
##  [769] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
##  [775] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
##  [781] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
##  [787] 1.986477110 1.986477110 1.986477110 1.986477110 0.061692868 0.061692868
##  [793] 0.061692868 0.061692868 0.061692868 0.074960750 0.074960750 0.088612338
##  [799] 0.088612338 0.088612338 0.088612338 0.088612338 0.088612338 0.117433167
##  [805] 0.132798388 0.132798388 0.132798388 0.132798388 0.148612854 0.148612854
##  [811] 0.148612854 0.148612854 0.148612854 0.165103160 0.165103160 0.165103160
##  [817] 0.181895734 0.199245046 0.199245046 0.199245046 0.199245046 0.199245046
##  [823] 0.217183883 0.217183883 0.217183883 0.217183883 0.235394131 0.255492719
##  [829] 0.276175004 0.276175004 0.276175004 0.297464973 0.297464973 0.297464973
##  [835] 0.345150450 0.345150450 0.369758625 0.369758625 0.369758625 0.395949724
##  [841] 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724
##  [847] 0.395949724 0.395949724 0.395949724 0.422705551 0.422705551 0.451879735
##  [853] 0.451879735 0.451879735 0.451879735 0.451879735 0.483048139 0.483048139
##  [859] 0.483048139 0.483048139 0.515749199 0.515749199 0.515749199 0.515749199
##  [865] 0.515749199 0.515749199 0.515749199 0.550363970 0.550363970 0.550363970
##  [871] 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
##  [877] 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
##  [883] 0.550363970 0.550363970 0.630418634 0.630418634 0.630418634 0.630418634
##  [889] 0.630418634 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [895] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [901] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [907] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [913] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [919] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [925] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [931] 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
##  [937] 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
##  [943] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924
##  [949] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924
##  [955] 0.775050924 0.775050924 0.837188229 0.837188229 0.837188229 0.837188229
##  [961] 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229
##  [967] 0.837188229 0.837188229 0.922115104 0.922115104 0.922115104 0.922115104
##  [973] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [979] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [985] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [991] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [997] 0.922115104 0.922115104 0.922115104 0.922115104 1.025522965 1.025522965
## [1003] 1.025522965 1.025522965 1.137150675 1.137150675 1.137150675 1.137150675
## [1009] 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675
## [1015] 1.137150675 1.137150675 1.269696765 1.269696765 1.269696765 1.269696765
## [1021] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1027] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1033] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1039] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1045] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1051] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1057] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1063] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1069] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1075] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.434236318
## [1081] 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318
## [1087] 1.640072935 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1093] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1099] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1105] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1111] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1117] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1123] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.781716629
## [1129] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1135] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1141] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1147] 2.781716629
##    [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
##    [7] 0.000000000 0.002422910 0.002422910 0.002422910 0.002422910 0.007342705
##   [13] 0.007342705 0.009818343 0.012299473 0.012299473 0.012299473 0.017311969
##   [19] 0.017311969 0.019844167 0.022381833 0.022381833 0.022381833 0.024936882
##   [25] 0.027499712 0.030075870 0.030075870 0.030075870 0.032669946 0.035270579
##   [31] 0.037874232 0.037874232 0.037874232 0.037874232 0.040499620 0.043138885
##   [37] 0.043138885 0.045803991 0.045803991 0.045803991 0.045803991 0.045803991
##   [43] 0.051234459 0.051234459 0.053959659 0.053959659 0.053959659 0.056709190
##   [49] 0.056709190 0.056709190 0.059476974 0.062246981 0.062246981 0.067836269
##   [55] 0.067836269 0.070640906 0.070640906 0.070640906 0.073473549 0.076333898
##   [61] 0.076333898 0.076333898 0.076333898 0.079228569 0.082130814 0.082130814
##   [67] 0.085048413 0.087975882 0.087975882 0.090949553 0.096957268 0.099985028
##   [73] 0.103028855 0.109139889 0.109139889 0.109139889 0.112228053 0.112228053
##   [79] 0.118449962 0.118449962 0.121585833 0.121585833 0.124740883 0.127910783
##   [85] 0.131086222 0.131086222 0.134295773 0.137512089 0.140741665 0.143980878
##   [91] 0.147245615 0.150524681 0.153811349 0.157132578 0.160467711 0.167198221
##   [97] 0.167198221 0.170591826 0.170591826 0.170591826 0.170591826 0.170591826
##  [103] 0.174035043 0.174035043 0.174035043 0.174035043 0.174035043 0.174035043
##  [109] 0.174035043 0.174035043 0.177557058 0.177557058 0.181106991 0.181106991
##  [115] 0.184688635 0.188276315 0.188276315 0.191893457 0.191893457 0.195545936
##  [121] 0.195545936 0.199220479 0.202918422 0.202918422 0.206640238 0.210388559
##  [127] 0.214153663 0.217940133 0.217940133 0.217940133 0.221759485 0.225583152
##  [133] 0.225583152 0.233296347 0.237181790 0.241074706 0.244987967 0.248918660
##  [139] 0.252879175 0.256857068 0.260841555 0.260841555 0.264850422 0.268867250
##  [145] 0.272931469 0.272931469 0.277008479 0.277008479 0.281097733 0.285209776
##  [151] 0.289338169 0.293505342 0.293505342 0.293505342 0.293505342 0.293505342
##  [157] 0.293505342 0.297817238 0.297817238 0.302187933 0.306586088 0.306586088
##  [163] 0.311044871 0.315508708 0.319986281 0.319986281 0.319986281 0.324531278
##  [169] 0.324531278 0.329130385 0.333778961 0.333778961 0.338472789 0.343182554
##  [175] 0.343182554 0.343182554 0.347921280 0.352686126 0.357475797 0.362277850
##  [181] 0.362277850 0.367093534 0.371936042 0.376814323 0.381700710 0.386607441
##  [187] 0.391531514 0.396487204 0.396487204 0.396487204 0.401486385 0.406493245
##  [193] 0.411533108 0.411533108 0.411533108 0.416642840 0.421783672 0.421783672
##  [199] 0.427027876 0.432287963 0.437574903 0.437574903 0.442910860 0.442910860
##  [205] 0.453714946 0.459186349 0.459186349 0.464725859 0.470285325 0.475873942
##  [211] 0.481490529 0.481490529 0.481490529 0.481490529 0.487209959 0.492988895
##  [217] 0.498820289 0.498820289 0.504785833 0.504785833 0.510816417 0.510816417
##  [223] 0.510816417 0.516962041 0.523119100 0.523119100 0.529304753 0.535538855
##  [229] 0.541833971 0.541833971 0.541833971 0.548306163 0.548306163 0.554827461
##  [235] 0.561383167 0.567992686 0.574625244 0.581315779 0.581315779 0.588074133
##  [241] 0.594886993 0.601782986 0.608740305 0.615751400 0.622848256 0.629999870
##  [247] 0.637258708 0.644540098 0.644540098 0.644540098 0.644540098 0.652030508
##  [253] 0.659604587 0.659604587 0.667301465 0.667301465 0.675122748 0.683036373
##  [259] 0.691009874 0.699068259 0.707188520 0.715354708 0.723609767 0.731903802
##  [265] 0.740285562 0.748838785 0.748838785 0.757639804 0.757639804 0.757639804
##  [271] 0.757639804 0.757639804 0.757639804 0.757639804 0.757639804 0.766967351
##  [277] 0.776382646 0.785932909 0.785932909 0.785932909 0.795662298 0.805630556
##  [283] 0.805630556 0.805630556 0.815715614 0.825948817 0.836324102 0.836324102
##  [289] 0.836324102 0.847059561 0.847059561 0.857909699 0.868808764 0.879793185
##  [295] 0.879793185 0.879793185 0.891329087 0.891329087 0.891329087 0.903058135
##  [301] 0.903058135 0.914948511 0.926971810 0.926971810 0.939203300 0.939203300
##  [307] 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300
##  [313] 0.939203300 0.939203300 0.939203300 0.952862217 0.966757789 0.966757789
##  [319] 0.966757789 0.966757789 0.981074755 0.995864031 0.995864031 1.010953971
##  [325] 1.010953971 1.010953971 1.026685796 1.042530764 1.042530764 1.042530764
##  [331] 1.042530764 1.042530764 1.059580899 1.076967045 1.076967045 1.076967045
##  [337] 1.076967045 1.076967045 1.076967045 1.076967045 1.076967045 1.076967045
##  [343] 1.116142086 1.116142086 1.136740519 1.136740519 1.158208569 1.180687663
##  [349] 1.180687663 1.204215922 1.204215922 1.204215922 1.229160124 1.254409050
##  [355] 1.280031424 1.280031424 1.306648637 1.333899572 1.333899572 0.000000000
##  [361] 0.000000000 0.000000000 0.000000000 0.000000000 0.012819623 0.012819623
##  [367] 0.026034998 0.039566407 0.039566407 0.039566407 0.053402500 0.053402500
##  [373] 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500
##  [379] 0.068122762 0.068122762 0.068122762 0.068122762 0.068122762 0.083749980
##  [385] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
##  [391] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
##  [397] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.100386010
##  [403] 0.100386010 0.100386010 0.100386010 0.100386010 0.100386010 0.117547315
##  [409] 0.134850559 0.134850559 0.134850559 0.134850559 0.134850559 0.171646666
##  [415] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
##  [421] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
##  [427] 0.171646666 0.171646666 0.171646666 0.191372504 0.191372504 0.191372504
##  [433] 0.191372504 0.191372504 0.191372504 0.211774184 0.211774184 0.211774184
##  [439] 0.211774184 0.211774184 0.211774184 0.211774184 0.232795610 0.232795610
##  [445] 0.232795610 0.232795610 0.232795610 0.232795610 0.276919329 0.276919329
##  [451] 0.276919329 0.276919329 0.300661572 0.300661572 0.300661572 0.300661572
##  [457] 0.300661572 0.300661572 0.300661572 0.324925481 0.324925481 0.324925481
##  [463] 0.324925481 0.324925481 0.350680784 0.350680784 0.350680784 0.350680784
##  [469] 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784
##  [475] 0.350680784 0.350680784 0.378294064 0.378294064 0.378294064 0.378294064
##  [481] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
##  [487] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
##  [493] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
##  [499] 0.406697290 0.406697290 0.406697290 0.436020806 0.436020806 0.436020806
##  [505] 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806
##  [511] 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806 0.466837198
##  [517] 0.466837198 0.466837198 0.466837198 0.466837198 0.466837198 0.466837198
##  [523] 0.466837198 0.466837198 0.466837198 0.499883881 0.499883881 0.499883881
##  [529] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
##  [535] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
##  [541] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.569580748
##  [547] 0.569580748 0.569580748 0.569580748 0.569580748 0.607417459 0.647833203
##  [553] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
##  [559] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
##  [565] 0.647833203 0.647833203 0.647833203 0.698345010 0.698345010 0.698345010
##  [571] 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010
##  [577] 0.698345010 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
##  [583] 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
##  [589] 0.756226472 0.756226472 0.756226472 0.818519693 0.818519693 0.818519693
##  [595] 0.818519693 0.889304890 0.889304890 0.889304890 0.889304890 0.889304890
##  [601] 0.889304890 0.968545149 0.968545149 0.968545149 1.060605655 1.060605655
##  [607] 1.060605655 1.060605655 1.060605655 1.060605655 1.060605655 1.162895503
##  [613] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
##  [619] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
##  [625] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
##  [631] 1.162895503 1.162895503 1.162895503 1.288056766 1.288056766 1.288056766
##  [637] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [643] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [649] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [655] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [661] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [667] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
##  [673] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [679] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [685] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [691] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [697] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
##  [703] 1.472408461 1.472408461 1.700899877 1.700899877 1.700899877 1.700899877
##  [709] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877
##  [715] 1.700899877 1.700899877 1.700899877 1.700899877 0.023953755 0.036281742
##  [721] 0.048854856 0.061692868 0.061692868 0.061692868 0.061692868 0.061692868
##  [727] 0.074960750 0.074960750 0.088612338 0.088612338 0.088612338 0.088612338
##  [733] 0.088612338 0.088612338 0.102901922 0.117433167 0.132798388 0.132798388
##  [739] 0.132798388 0.132798388 0.148612854 0.148612854 0.148612854 0.148612854
##  [745] 0.148612854 0.165103160 0.165103160 0.165103160 0.165103160 0.181895734
##  [751] 0.181895734 0.199245046 0.199245046 0.199245046 0.199245046 0.199245046
##  [757] 0.217183883 0.217183883 0.217183883 0.217183883 0.235394131 0.255492719
##  [763] 0.255492719 0.276175004 0.276175004 0.276175004 0.276175004 0.297464973
##  [769] 0.297464973 0.297464973 0.320944085 0.345150450 0.369758625 0.369758625
##  [775] 0.369758625 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724
##  [781] 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724
##  [787] 0.422705551 0.422705551 0.451879735 0.451879735 0.451879735 0.451879735
##  [793] 0.451879735 0.483048139 0.483048139 0.483048139 0.483048139 0.515749199
##  [799] 0.515749199 0.515749199 0.515749199 0.515749199 0.515749199 0.515749199
##  [805] 0.515749199 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
##  [811] 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
##  [817] 0.550363970 0.550363970 0.589208062 0.630418634 0.630418634 0.630418634
##  [823] 0.630418634 0.630418634 0.673825328 0.673825328 0.673825328 0.673825328
##  [829] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [835] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [841] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [847] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [853] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [859] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
##  [865] 0.673825328 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
##  [871] 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
##  [877] 0.720528533 0.720528533 0.720528533 0.775050924 0.775050924 0.775050924
##  [883] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924
##  [889] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.837188229
##  [895] 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229
##  [901] 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229
##  [907] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [913] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [919] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [925] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
##  [931] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 1.025522965
##  [937] 1.025522965 1.025522965 1.025522965 1.137150675 1.137150675 1.137150675
##  [943] 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675
##  [949] 1.137150675 1.137150675 1.269696765 1.269696765 1.269696765 1.269696765
##  [955] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
##  [961] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
##  [967] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
##  [973] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
##  [979] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
##  [985] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
##  [991] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
##  [997] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1003] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1009] 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318
## [1015] 1.434236318 1.434236318 1.640072935 1.640072935 2.013928967 2.013928967
## [1021] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1027] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1033] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1039] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1045] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1051] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1057] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1063] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1069] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1075] 2.781716629 2.781716629 2.781716629
##   [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.002422910
##   [7] 0.002422910 0.002422910 0.002422910 0.002422910 0.009818343 0.012299473
##  [13] 0.012299473 0.017311969 0.017311969 0.019844167 0.022381833 0.022381833
##  [19] 0.024936882 0.030075870 0.037874232 0.037874232 0.045803991 0.051234459
##  [25] 0.051234459 0.053959659 0.056709190 0.056709190 0.062246981 0.067836269
##  [31] 0.067836269 0.067836269 0.070640906 0.087975882 0.096957268 0.109139889
##  [37] 0.109139889 0.118449962 0.127910783 0.137512089 0.140741665 0.140741665
##  [43] 0.153811349 0.170591826 0.170591826 0.170591826 0.174035043 0.174035043
##  [49] 0.174035043 0.188276315 0.199220479 0.260841555 0.277008479 0.293505342
##  [55] 0.293505342 0.306586088 0.319986281 0.324531278 0.343182554 0.343182554
##  [61] 0.362277850 0.386607441 0.396487204 0.406493245 0.406493245 0.411533108
##  [67] 0.464725859 0.481490529 0.481490529 0.498820289 0.498820289 0.510816417
##  [73] 0.535538855 0.548306163 0.548306163 0.567992686 0.601782986 0.622848256
##  [79] 0.652030508 0.757639804 0.757639804 0.847059561 0.891329087 0.903058135
##  [85] 0.903058135 0.939203300 0.939203300 1.076967045 1.136740519 1.333899572
##  [91] 0.000000000 0.000000000 0.000000000 0.000000000 0.012819623 0.026034998
##  [97] 0.039566407 0.039566407 0.039566407 0.053402500 0.053402500 0.053402500
## [103] 0.053402500 0.053402500 0.053402500 0.068122762 0.068122762 0.068122762
## [109] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.100386010
## [115] 0.100386010 0.100386010 0.117547315 0.134850559 0.134850559 0.134850559
## [121] 0.153076521 0.171646666 0.171646666 0.171646666 0.191372504 0.191372504
## [127] 0.211774184 0.211774184 0.232795610 0.232795610 0.254560972 0.276919329
## [133] 0.300661572 0.324925481 0.324925481 0.324925481 0.350680784 0.350680784
## [139] 0.350680784 0.378294064 0.406697290 0.436020806 0.436020806 0.436020806
## [145] 0.466837198 0.466837198 0.499883881 0.534298268 0.534298268 0.534298268
## [151] 0.534298268 0.569580748 0.569580748 0.607417459 0.647833203 0.647833203
## [157] 0.698345010 0.698345010 0.698345010 0.698345010 0.756226472 0.756226472
## [163] 0.756226472 0.818519693 0.818519693 0.889304890 0.968545149 1.060605655
## [169] 1.162895503 1.162895503 1.288056766 1.288056766 1.288056766 1.288056766
## [175] 1.288056766 1.472408461 1.472408461 1.472408461 1.700899877 1.700899877
## [181] 0.023953755 0.036281742 0.048854856 0.061692868 0.061692868 0.061692868
## [187] 0.074960750 0.074960750 0.088612338 0.088612338 0.088612338 0.088612338
## [193] 0.102901922 0.117433167 0.132798388 0.132798388 0.132798388 0.148612854
## [199] 0.148612854 0.165103160 0.181895734 0.199245046 0.217183883 0.235394131
## [205] 0.255492719 0.276175004 0.276175004 0.297464973 0.320944085 0.345150450
## [211] 0.345150450 0.369758625 0.395949724 0.422705551 0.451879735 0.483048139
## [217] 0.483048139 0.515749199 0.515749199 0.550363970 0.550363970 0.550363970
## [223] 0.550363970 0.550363970 0.589208062 0.630418634 0.630418634 0.673825328
## [229] 0.673825328 0.673825328 0.673825328 0.673825328 0.720528533 0.720528533
## [235] 0.720528533 0.775050924 0.775050924 0.775050924 0.775050924 0.837188229
## [241] 0.837188229 0.837188229 0.922115104 0.922115104 0.922115104 0.922115104
## [247] 0.922115104 0.922115104 0.922115104 1.025522965 1.025522965 1.137150675
## [253] 1.137150675 1.137150675 1.269696765 1.269696765 1.269696765 1.269696765
## [259] 1.269696765 1.269696765 1.269696765 1.434236318 1.640072935 2.013928967
## [265] 2.013928967 2.013928967 2.013928967 2.781716629 2.781716629 2.781716629
res3<-residuals(Cox.fit3,type="deviance")
hist(res3)
<<<<<<< HEAD

boxplot(res3~Cdata$Death_Censoring)

=======

boxplot(res3~Cdata$Death_Censoring)

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
# Looking at Age
min(Cdata[,7])# youngest person is 17
## [1] 17
max(Cdata[,7])# oldest person is 80
## [1] 80
Age.groups=cut(Cdata$Age,breaks=c(-Inf,quantile(Cdata$Age ,prob=c(0.33,0.66)),Inf)) #devides the Age into groups
Cox.fit4 <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ strata(Age.groups)+ Prothrombin, data = Cdata)
summary(Cox.fit4)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     strata(Age.groups) + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                         coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Treatment)1  0.026418  1.026770  0.119345  0.221   0.8248    
## factor(Sex)1        0.497157  1.644041  0.127598  3.896 9.77e-05 ***
## factor(Ascites)1    0.459189  1.582789  0.183471  2.503   0.0123 *  
## factor(Ascites)2    1.046229  2.846896  0.184316  5.676 1.38e-08 ***
## Prothrombin        -0.012052  0.988020  0.002876 -4.191 2.78e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                    exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1     1.027     0.9739    0.8126    1.2974
## factor(Sex)1           1.644     0.6083    1.2803    2.1112
## factor(Ascites)1       1.583     0.6318    1.1047    2.2677
## factor(Ascites)2       2.847     0.3513    1.9837    4.0857
## Prothrombin            0.988     1.0121    0.9825    0.9936
## 
## Concordance= 0.641  (se = 0.019 )
## Likelihood ratio test= 73.49  on 5 df,   p=2e-14
## Wald test            = 78.88  on 5 df,   p=1e-15
## Score (logrank) test = 84.54  on 5 df,   p=<2e-16
source("andersen.R")
test4<-survfit(Cox.fit4)
test4
## Call: survfit(formula = Cox.fit4)
## 
##               n events median 0.95LCL 0.95UCL
## (-Inf,55.7] 161     68   3078    2467    3611
## (55.7,64]   169    106   1559    1263    2057
## (64, Inf]   158    118    852     572    1169
andersen.plot(test4) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD

=======

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##   [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.005528288
##   [7] 0.005528288 0.005528288 0.005528288 0.005528288 0.005528288 0.005528288
##  [13] 0.005528288 0.011403634 0.017304927 0.017304927 0.017304927 0.017304927
##  [19] 0.017304927 0.023338937 0.029409752 0.029409752 0.029409752 0.035588296
##  [25] 0.035588296 0.035588296 0.041919831 0.041919831 0.048303041 0.048303041
##  [31] 0.048303041 0.054918882 0.061553184 0.068375151 0.068375151 0.068375151
##  [37] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
##  [43] 0.068375151 0.075365321 0.075365321 0.075365321 0.082455877 0.082455877
##  [49] 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877
##  [55] 0.082455877 0.082455877 0.089674197 0.096968283 0.096968283 0.104306718
##  [61] 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718 0.111695071
##  [67] 0.111695071 0.111695071 0.119142818 0.126825522 0.126825522 0.126825522
##  [73] 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522
##  [79] 0.126825522 0.134984552 0.134984552 0.143227902 0.143227902 0.143227902
##  [85] 0.151634521 0.151634521 0.151634521 0.151634521 0.151634521 0.151634521
##  [91] 0.151634521 0.160282070 0.160282070 0.168962279 0.168962279 0.177720498
##  [97] 0.177720498 0.177720498 0.186541816 0.195438275 0.195438275 0.195438275
## [103] 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [109] 0.195438275 0.204585450 0.214054867 0.214054867 0.214054867 0.214054867
## [115] 0.223603523 0.223603523 0.233381313 0.243333679 0.243333679 0.253381794
## [121] 0.253381794 0.263515034 0.263515034 0.263515034 0.273710658 0.273710658
## [127] 0.273710658 0.273710658 0.283960762 0.283960762 0.283960762 0.283960762
## [133] 0.283960762 0.283960762 0.283960762 0.294431829 0.294431829 0.305000280
## [139] 0.305000280 0.305000280 0.305000280 0.305000280 0.305000280 0.305000280
## [145] 0.305000280 0.315826690 0.315826690 0.315826690 0.326983499 0.326983499
## [151] 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499 0.338812149
## [157] 0.338812149 0.338812149 0.338812149 0.350978763 0.363224224 0.363224224
## [163] 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224
## [169] 0.363224224 0.363224224 0.363224224 0.376029567 0.376029567 0.376029567
## [175] 0.376029567 0.376029567 0.389106667 0.389106667 0.402380512 0.416010968
## [181] 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813
## [187] 0.429762813 0.444009105 0.444009105 0.444009105 0.458561726 0.473262284
## [193] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [199] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.489049109
## [205] 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109
## [211] 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109 0.505600045
## [217] 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045
## [223] 0.505600045 0.523987973 0.523987973 0.543062495 0.543062495 0.543062495
## [229] 0.543062495 0.543062495 0.543062495 0.543062495 0.563339269 0.563339269
## [235] 0.563339269 0.563339269 0.586028574 0.586028574 0.586028574 0.586028574
## [241] 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574
## [247] 0.586028574 0.586028574 0.586028574 0.612099303 0.612099303 0.612099303
## [253] 0.612099303 0.612099303 0.639620829 0.639620829 0.639620829 0.639620829
## [259] 0.670502874 0.670502874 0.670502874 0.702589472 0.702589472 0.737090896
## [265] 0.737090896 0.737090896 0.737090896 0.777173206 0.818773352 0.818773352
## [271] 0.818773352 0.861528837 0.906277419 0.906277419 0.906277419 0.906277419
## [277] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [283] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [289] 0.959570305 0.959570305 0.959570305 1.090085142 1.090085142 1.090085142
## [295] 1.090085142 1.090085142 1.090085142 1.090085142 1.090085142 1.090085142
## [301] 1.090085142 1.090085142 1.090085142 0.005350627 0.010822944 0.010822944
## [307] 0.010822944 0.016454369 0.016454369 0.016454369 0.022164713 0.022164713
## [313] 0.027889662 0.027889662 0.027889662 0.033763774 0.033763774 0.033763774
## [319] 0.033763774 0.033763774 0.033763774 0.039791244 0.039791244 0.039791244
## [325] 0.039791244 0.045954956 0.045954956 0.045954956 0.045954956 0.045954956
## [331] 0.052253850 0.052253850 0.052253850 0.058725586 0.058725586 0.058725586
## [337] 0.065226295 0.065226295 0.071800183 0.071800183 0.071800183 0.071800183
## [343] 0.078528422 0.078528422 0.085296520 0.085296520 0.085296520 0.092181968
## [349] 0.099100984 0.099100984 0.106076035 0.106076035 0.113136184 0.120399425
## [355] 0.120399425 0.127705587 0.127705587 0.135076428 0.142522708 0.142522708
## [361] 0.142522708 0.149993842 0.149993842 0.157531961 0.165150610 0.165150610
## [367] 0.172846201 0.180620008 0.188445918 0.188445918 0.196368910 0.196368910
## [373] 0.196368910 0.196368910 0.196368910 0.196368910 0.204695686 0.204695686
## [379] 0.204695686 0.204695686 0.204695686 0.213256293 0.213256293 0.213256293
## [385] 0.213256293 0.221921101 0.221921101 0.221921101 0.230668576 0.230668576
## [391] 0.239491574 0.248425299 0.248425299 0.248425299 0.248425299 0.248425299
## [397] 0.257512175 0.266684841 0.266684841 0.266684841 0.275887140 0.275887140
## [403] 0.275887140 0.275887140 0.285301563 0.294762418 0.304319540 0.304319540
## [409] 0.304319540 0.304319540 0.304319540 0.314131316 0.314131316 0.314131316
## [415] 0.324057279 0.334041522 0.334041522 0.334041522 0.344145476 0.344145476
## [421] 0.344145476 0.344145476 0.344145476 0.354489364 0.354489364 0.364955963
## [427] 0.375468287 0.375468287 0.386067530 0.396810622 0.396810622 0.396810622
## [433] 0.418872475 0.430217353 0.430217353 0.430217353 0.441771494 0.453445406
## [439] 0.453445406 0.465319253 0.465319253 0.465319253 0.477384731 0.477384731
## [445] 0.489634257 0.502028814 0.502028814 0.514632492 0.514632492 0.527340450
## [451] 0.540164985 0.540164985 0.540164985 0.540164985 0.540164985 0.553304597
## [457] 0.566657491 0.580311979 0.580311979 0.580311979 0.594305960 0.594305960
## [463] 0.594305960 0.594305960 0.594305960 0.609012152 0.623904225 0.623904225
## [469] 0.639143157 0.639143157 0.654807088 0.670718435 0.687301167 0.704071054
## [475] 0.704071054 0.721157987 0.721157987 0.738561528 0.756346779 0.756346779
## [481] 0.774488197 0.774488197 0.774488197 0.774488197 0.774488197 0.793145557
## [487] 0.812226345 0.812226345 0.812226345 0.831644405 0.831644405 0.851361126
## [493] 0.871448014 0.871448014 0.871448014 0.891810367 0.891810367 0.913268483
## [499] 0.913268483 0.913268483 0.913268483 0.913268483 0.913268483 0.913268483
## [505] 0.913268483 0.936798339 0.936798339 0.936798339 0.936798339 0.961841394
## [511] 0.987712859 0.987712859 0.987712859 1.013960311 1.040953046 1.040953046
## [517] 1.040953046 1.069941379 1.069941379 1.099543534 1.129531930 1.129531930
## [523] 1.160547040 1.192692166 1.192692166 1.228912605 1.228912605 1.228912605
## [529] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [535] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [541] 1.228912605 1.228912605 1.228912605 1.270816176 1.270816176 1.314046197
## [547] 1.314046197 1.359029219 1.404934747 1.404934747 1.404934747 1.404934747
## [553] 1.404934747 1.458863406 1.458863406 1.458863406 1.458863406 1.458863406
## [559] 1.458863406 1.458863406 1.458863406 1.520684154 1.520684154 1.590254427
## [565] 1.590254427 1.590254427 1.590254427 1.590254427 1.667260347 1.667260347
## [571] 1.667260347 1.667260347 1.747291122 1.747291122 1.747291122 1.747291122
## [577] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [583] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [589] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.926625374
## [595] 1.926625374 1.926625374 1.926625374 1.926625374 1.926625374 2.281783691
## [601] 2.281783691 2.281783691 2.281783691 2.281783691 2.281783691 2.281783691
## [607] 0.009633169 0.014551755 0.019546218 0.019546218 0.024599344 0.029678509
## [613] 0.029678509 0.045413686 0.045413686 0.045413686 0.050750288 0.050750288
## [619] 0.050750288 0.061572580 0.072742518 0.078415489 0.078415489 0.107873484
## [625] 0.107873484 0.107873484 0.107873484 0.126231259 0.126231259 0.138970374
## [631] 0.138970374 0.138970374 0.158792804 0.158792804 0.165618521 0.165618521
## [637] 0.172484282 0.186450350 0.193667471 0.200947081 0.230731258 0.230731258
## [643] 0.246241191 0.246241191 0.246241191 0.246241191 0.254091780 0.254091780
## [649] 0.269933185 0.269933185 0.295740183 0.295740183 0.295740183 0.304607658
## [655] 0.304607658 0.304607658 0.313541365 0.313541365 0.322590148 0.341155661
## [661] 0.341155661 0.350635123 0.350635123 0.350635123 0.360361864 0.360361864
## [667] 0.360361864 0.391154668 0.401655837 0.401655837 0.401655837 0.412274302
## [673] 0.412274302 0.433912349 0.433912349 0.433912349 0.433912349 0.433912349
## [679] 0.445124076 0.445124076 0.445124076 0.445124076 0.445124076 0.456660219
## [685] 0.456660219 0.456660219 0.468424656 0.468424656 0.480317909 0.480317909
## [691] 0.492387762 0.492387762 0.492387762 0.492387762 0.504593246 0.516909136
## [697] 0.516909136 0.516909136 0.529347066 0.529347066 0.529347066 0.580175913
## [703] 0.593186674 0.593186674 0.593186674 0.606344153 0.606344153 0.606344153
## [709] 0.619675763 0.619675763 0.647161102 0.647161102 0.661946165 0.677209373
## [715] 0.677209373 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700
## [721] 0.692791700 0.692791700 0.692791700 0.708979609 0.708979609 0.708979609
## [727] 0.742673869 0.760665334 0.760665334 0.760665334 0.779916835 0.839297067
## [733] 0.839297067 0.859999536 0.859999536 0.859999536 0.881788673 0.881788673
## [739] 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673
## [745] 0.881788673 0.881788673 0.881788673 0.905686340 0.905686340 0.905686340
## [751] 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340
## [757] 0.905686340 0.905686340 0.930593182 0.956024835 0.956024835 0.956024835
## [763] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [769] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [775] 0.982860461 1.011138349 1.011138349 1.011138349 1.072147856 1.072147856
## [781] 1.104113664 1.104113664 1.136775819 1.171345121 1.171345121 1.207040278
## [787] 1.207040278 1.207040278 1.207040278 1.246907275 1.246907275 1.287930715
## [793] 1.330411123 1.330411123 1.330411123 1.330411123 1.374604413 1.374604413
## [799] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943
## [805] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.521709021
## [811] 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021
## [817] 1.574890670 1.574890670 1.574890670 1.574890670 1.574890670 1.574890670
## [823] 1.574890670 1.685791093 1.685791093 1.750311642 1.750311642 1.750311642
## [829] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [835] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [841] 1.750311642 1.750311642 1.832400858 1.832400858 1.832400858 1.832400858
## [847] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [853] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [859] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [865] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 2.029441753
## [871] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [877] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [883] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.194031226
## [889] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226
## [895] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.530841239
## [901] 2.925483659 2.925483659 2.925483659
##   [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
##   [7] 0.000000000 0.005528288 0.005528288 0.005528288 0.005528288 0.005528288
##  [13] 0.005528288 0.005528288 0.005528288 0.011403634 0.011403634 0.017304927
##  [19] 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927
##  [25] 0.017304927 0.023338937 0.029409752 0.029409752 0.029409752 0.029409752
##  [31] 0.029409752 0.029409752 0.029409752 0.029409752 0.035588296 0.035588296
##  [37] 0.041919831 0.041919831 0.048303041 0.048303041 0.048303041 0.048303041
##  [43] 0.054918882 0.054918882 0.061553184 0.068375151 0.068375151 0.068375151
##  [49] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
##  [55] 0.068375151 0.068375151 0.075365321 0.075365321 0.082455877 0.082455877
##  [61] 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877
##  [67] 0.082455877 0.089674197 0.089674197 0.096968283 0.096968283 0.104306718
##  [73] 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718 0.111695071
##  [79] 0.111695071 0.111695071 0.111695071 0.111695071 0.111695071 0.119142818
##  [85] 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522
##  [91] 0.126825522 0.126825522 0.126825522 0.134984552 0.134984552 0.134984552
##  [97] 0.143227902 0.143227902 0.143227902 0.143227902 0.151634521 0.151634521
## [103] 0.151634521 0.151634521 0.151634521 0.160282070 0.160282070 0.168962279
## [109] 0.168962279 0.168962279 0.168962279 0.168962279 0.177720498 0.177720498
## [115] 0.177720498 0.186541816 0.186541816 0.195438275 0.195438275 0.195438275
## [121] 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [127] 0.195438275 0.195438275 0.204585450 0.214054867 0.214054867 0.223603523
## [133] 0.223603523 0.223603523 0.233381313 0.233381313 0.243333679 0.253381794
## [139] 0.253381794 0.253381794 0.253381794 0.253381794 0.263515034 0.263515034
## [145] 0.263515034 0.273710658 0.273710658 0.273710658 0.273710658 0.273710658
## [151] 0.283960762 0.283960762 0.283960762 0.283960762 0.283960762 0.294431829
## [157] 0.305000280 0.305000280 0.305000280 0.305000280 0.315826690 0.326983499
## [163] 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499 0.338812149
## [169] 0.338812149 0.338812149 0.338812149 0.350978763 0.363224224 0.363224224
## [175] 0.363224224 0.363224224 0.363224224 0.376029567 0.376029567 0.376029567
## [181] 0.376029567 0.376029567 0.376029567 0.389106667 0.389106667 0.402380512
## [187] 0.402380512 0.416010968 0.429762813 0.429762813 0.429762813 0.429762813
## [193] 0.429762813 0.429762813 0.429762813 0.429762813 0.444009105 0.444009105
## [199] 0.458561726 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [205] 0.473262284 0.473262284 0.473262284 0.489049109 0.489049109 0.489049109
## [211] 0.489049109 0.489049109 0.489049109 0.505600045 0.505600045 0.505600045
## [217] 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045 0.523987973
## [223] 0.523987973 0.543062495 0.543062495 0.543062495 0.543062495 0.543062495
## [229] 0.543062495 0.563339269 0.563339269 0.563339269 0.563339269 0.586028574
## [235] 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574
## [241] 0.586028574 0.586028574 0.612099303 0.612099303 0.612099303 0.612099303
## [247] 0.639620829 0.639620829 0.639620829 0.670502874 0.702589472 0.702589472
## [253] 0.702589472 0.737090896 0.737090896 0.737090896 0.737090896 0.777173206
## [259] 0.818773352 0.861528837 0.906277419 0.906277419 0.906277419 0.959570305
## [265] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [271] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [277] 0.959570305 1.090085142 1.090085142 1.090085142 1.090085142 1.090085142
## [283] 1.090085142 1.090085142 0.000000000 0.005350627 0.010822944 0.010822944
## [289] 0.010822944 0.010822944 0.016454369 0.016454369 0.016454369 0.016454369
## [295] 0.016454369 0.022164713 0.022164713 0.027889662 0.027889662 0.033763774
## [301] 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774
## [307] 0.033763774 0.033763774 0.033763774 0.039791244 0.039791244 0.039791244
## [313] 0.039791244 0.039791244 0.039791244 0.045954956 0.045954956 0.045954956
## [319] 0.045954956 0.045954956 0.045954956 0.052253850 0.052253850 0.052253850
## [325] 0.058725586 0.058725586 0.058725586 0.058725586 0.058725586 0.065226295
## [331] 0.065226295 0.065226295 0.065226295 0.071800183 0.071800183 0.071800183
## [337] 0.078528422 0.078528422 0.085296520 0.085296520 0.085296520 0.092181968
## [343] 0.099100984 0.106076035 0.113136184 0.120399425 0.127705587 0.127705587
## [349] 0.127705587 0.127705587 0.142522708 0.142522708 0.142522708 0.142522708
## [355] 0.149993842 0.149993842 0.157531961 0.157531961 0.157531961 0.165150610
## [361] 0.180620008 0.188445918 0.188445918 0.188445918 0.188445918 0.188445918
## [367] 0.188445918 0.196368910 0.196368910 0.196368910 0.196368910 0.196368910
## [373] 0.204695686 0.204695686 0.204695686 0.204695686 0.213256293 0.213256293
## [379] 0.213256293 0.213256293 0.213256293 0.213256293 0.213256293 0.221921101
## [385] 0.221921101 0.230668576 0.239491574 0.248425299 0.248425299 0.248425299
## [391] 0.248425299 0.257512175 0.266684841 0.266684841 0.266684841 0.266684841
## [397] 0.266684841 0.266684841 0.266684841 0.275887140 0.275887140 0.275887140
## [403] 0.275887140 0.285301563 0.304319540 0.304319540 0.304319540 0.304319540
## [409] 0.304319540 0.304319540 0.304319540 0.314131316 0.314131316 0.314131316
## [415] 0.334041522 0.334041522 0.334041522 0.344145476 0.344145476 0.344145476
## [421] 0.344145476 0.344145476 0.344145476 0.344145476 0.344145476 0.354489364
## [427] 0.354489364 0.364955963 0.375468287 0.375468287 0.375468287 0.375468287
## [433] 0.375468287 0.396810622 0.396810622 0.418872475 0.430217353 0.430217353
## [439] 0.453445406 0.453445406 0.465319253 0.477384731 0.489634257 0.502028814
## [445] 0.514632492 0.540164985 0.540164985 0.540164985 0.540164985 0.540164985
## [451] 0.553304597 0.580311979 0.580311979 0.580311979 0.594305960 0.594305960
## [457] 0.594305960 0.594305960 0.623904225 0.639143157 0.670718435 0.704071054
## [463] 0.704071054 0.704071054 0.721157987 0.721157987 0.756346779 0.756346779
## [469] 0.756346779 0.774488197 0.774488197 0.774488197 0.774488197 0.774488197
## [475] 0.774488197 0.793145557 0.812226345 0.812226345 0.812226345 0.812226345
## [481] 0.831644405 0.871448014 0.871448014 0.871448014 0.871448014 0.913268483
## [487] 0.913268483 0.913268483 0.936798339 0.936798339 0.936798339 0.936798339
## [493] 0.936798339 0.961841394 0.987712859 0.987712859 0.987712859 1.069941379
## [499] 1.099543534 1.099543534 1.099543534 1.099543534 1.129531930 1.129531930
## [505] 1.192692166 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [511] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [517] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.270816176
## [523] 1.314046197 1.314046197 1.404934747 1.404934747 1.404934747 1.458863406
## [529] 1.458863406 1.458863406 1.458863406 1.458863406 1.458863406 1.520684154
## [535] 1.590254427 1.590254427 1.590254427 1.590254427 1.590254427 1.667260347
## [541] 1.667260347 1.667260347 1.667260347 1.747291122 1.747291122 1.747291122
## [547] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [553] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [559] 1.747291122 1.747291122 1.747291122 1.926625374 1.926625374 1.926625374
## [565] 1.926625374 1.926625374 1.926625374 2.281783691 0.009633169 0.014551755
## [571] 0.014551755 0.019546218 0.019546218 0.019546218 0.024599344 0.029678509
## [577] 0.029678509 0.034864290 0.040107303 0.045413686 0.045413686 0.050750288
## [583] 0.050750288 0.061572580 0.067117647 0.072742518 0.078415489 0.078415489
## [589] 0.084121378 0.089991681 0.095877420 0.101853535 0.107873484 0.107873484
## [595] 0.107873484 0.113924256 0.119989528 0.119989528 0.126231259 0.132564626
## [601] 0.132564626 0.138970374 0.138970374 0.138970374 0.158792804 0.165618521
## [607] 0.165618521 0.172484282 0.172484282 0.186450350 0.186450350 0.193667471
## [613] 0.193667471 0.200947081 0.208250442 0.223100168 0.230731258 0.238451311
## [619] 0.246241191 0.246241191 0.254091780 0.254091780 0.261998834 0.269933185
## [625] 0.269933185 0.295740183 0.295740183 0.304607658 0.313541365 0.313541365
## [631] 0.322590148 0.331761267 0.341155661 0.341155661 0.350635123 0.350635123
## [637] 0.350635123 0.350635123 0.360361864 0.360361864 0.370472540 0.380732678
## [643] 0.391154668 0.401655837 0.412274302 0.412274302 0.412274302 0.422970130
## [649] 0.422970130 0.422970130 0.433912349 0.433912349 0.433912349 0.433912349
## [655] 0.445124076 0.445124076 0.445124076 0.445124076 0.456660219 0.456660219
## [661] 0.456660219 0.456660219 0.468424656 0.468424656 0.468424656 0.480317909
## [667] 0.480317909 0.480317909 0.492387762 0.492387762 0.504593246 0.516909136
## [673] 0.516909136 0.516909136 0.529347066 0.529347066 0.541876644 0.554515396
## [679] 0.567272570 0.580175913 0.580175913 0.593186674 0.593186674 0.593186674
## [685] 0.606344153 0.606344153 0.606344153 0.619675763 0.633344103 0.647161102
## [691] 0.647161102 0.647161102 0.647161102 0.661946165 0.677209373 0.692791700
## [697] 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700
## [703] 0.692791700 0.708979609 0.708979609 0.708979609 0.725379668 0.725379668
## [709] 0.742673869 0.760665334 0.760665334 0.760665334 0.779916835 0.779916835
## [715] 0.799534319 0.819325249 0.839297067 0.859999536 0.859999536 0.881788673
## [721] 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673
## [727] 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340
## [733] 0.905686340 0.930593182 0.956024835 0.956024835 0.956024835 0.982860461
## [739] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [745] 1.011138349 1.011138349 1.040995345 1.072147856 1.072147856 1.104113664
## [751] 1.136775819 1.136775819 1.171345121 1.171345121 1.207040278 1.207040278
## [757] 1.207040278 1.207040278 1.207040278 1.246907275 1.246907275 1.287930715
## [763] 1.287930715 1.330411123 1.330411123 1.374604413 1.374604413 1.374604413
## [769] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.470281401
## [775] 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021
## [781] 1.574890670 1.574890670 1.629630555 1.685791093 1.685791093 1.685791093
## [787] 1.685791093 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [793] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [799] 1.750311642 1.750311642 1.832400858 1.832400858 1.832400858 1.832400858
## [805] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [811] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [817] 1.917615102 1.917615102 1.917615102 1.917615102 2.029441753 2.029441753
## [823] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [829] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [835] 2.029441753 2.029441753 2.029441753 2.194031226 2.194031226 2.194031226
## [841] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226
## [847] 2.194031226 2.530841239 2.925483659 2.925483659 2.925483659 2.925483659
##   [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
##   [7] 0.000000000 0.000000000 0.005528288 0.005528288 0.005528288 0.005528288
##  [13] 0.005528288 0.005528288 0.005528288 0.011403634 0.011403634 0.017304927
##  [19] 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927
##  [25] 0.017304927 0.017304927 0.017304927 0.029409752 0.029409752 0.029409752
##  [31] 0.029409752 0.029409752 0.029409752 0.029409752 0.029409752 0.035588296
##  [37] 0.035588296 0.035588296 0.041919831 0.041919831 0.041919831 0.048303041
##  [43] 0.048303041 0.048303041 0.048303041 0.054918882 0.068375151 0.068375151
##  [49] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
##  [55] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
##  [61] 0.068375151 0.075365321 0.075365321 0.075365321 0.082455877 0.082455877
##  [67] 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877
##  [73] 0.082455877 0.082455877 0.089674197 0.096968283 0.104306718 0.104306718
##  [79] 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718
##  [85] 0.104306718 0.111695071 0.111695071 0.111695071 0.111695071 0.111695071
##  [91] 0.111695071 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522
##  [97] 0.126825522 0.126825522 0.126825522 0.134984552 0.143227902 0.143227902
## [103] 0.143227902 0.151634521 0.151634521 0.151634521 0.151634521 0.151634521
## [109] 0.151634521 0.151634521 0.160282070 0.160282070 0.168962279 0.168962279
## [115] 0.168962279 0.168962279 0.168962279 0.168962279 0.177720498 0.177720498
## [121] 0.186541816 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [127] 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [133] 0.195438275 0.195438275 0.195438275 0.214054867 0.214054867 0.223603523
## [139] 0.223603523 0.223603523 0.233381313 0.243333679 0.253381794 0.253381794
## [145] 0.253381794 0.253381794 0.253381794 0.263515034 0.263515034 0.263515034
## [151] 0.263515034 0.273710658 0.273710658 0.273710658 0.273710658 0.273710658
## [157] 0.273710658 0.273710658 0.283960762 0.283960762 0.283960762 0.283960762
## [163] 0.283960762 0.283960762 0.283960762 0.283960762 0.294431829 0.305000280
## [169] 0.305000280 0.305000280 0.305000280 0.305000280 0.315826690 0.315826690
## [175] 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499
## [181] 0.326983499 0.338812149 0.338812149 0.338812149 0.363224224 0.363224224
## [187] 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224
## [193] 0.363224224 0.363224224 0.376029567 0.376029567 0.376029567 0.376029567
## [199] 0.376029567 0.376029567 0.376029567 0.389106667 0.389106667 0.402380512
## [205] 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813
## [211] 0.429762813 0.429762813 0.429762813 0.444009105 0.444009105 0.444009105
## [217] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [223] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [229] 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109
## [235] 0.489049109 0.489049109 0.489049109 0.505600045 0.505600045 0.505600045
## [241] 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045
## [247] 0.505600045 0.543062495 0.543062495 0.543062495 0.543062495 0.543062495
## [253] 0.543062495 0.563339269 0.563339269 0.586028574 0.586028574 0.586028574
## [259] 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574
## [265] 0.586028574 0.612099303 0.612099303 0.612099303 0.612099303 0.639620829
## [271] 0.670502874 0.670502874 0.670502874 0.702589472 0.737090896 0.737090896
## [277] 0.818773352 0.818773352 0.906277419 0.906277419 0.906277419 0.959570305
## [283] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [289] 0.959570305 0.959570305 1.090085142 1.090085142 1.090085142 1.090085142
## [295] 1.090085142 1.090085142 0.000000000 0.005350627 0.005350627 0.010822944
## [301] 0.010822944 0.010822944 0.010822944 0.016454369 0.016454369 0.016454369
## [307] 0.016454369 0.022164713 0.027889662 0.027889662 0.033763774 0.033763774
## [313] 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774
## [319] 0.033763774 0.033763774 0.033763774 0.033763774 0.039791244 0.039791244
## [325] 0.039791244 0.039791244 0.039791244 0.045954956 0.045954956 0.045954956
## [331] 0.045954956 0.045954956 0.045954956 0.045954956 0.045954956 0.052253850
## [337] 0.052253850 0.052253850 0.052253850 0.058725586 0.058725586 0.058725586
## [343] 0.065226295 0.065226295 0.065226295 0.065226295 0.065226295 0.071800183
## [349] 0.071800183 0.071800183 0.071800183 0.071800183 0.078528422 0.078528422
## [355] 0.085296520 0.085296520 0.085296520 0.085296520 0.092181968 0.099100984
## [361] 0.106076035 0.106076035 0.113136184 0.113136184 0.120399425 0.127705587
## [367] 0.127705587 0.127705587 0.135076428 0.142522708 0.142522708 0.149993842
## [373] 0.157531961 0.157531961 0.157531961 0.157531961 0.165150610 0.165150610
## [379] 0.172846201 0.180620008 0.180620008 0.188445918 0.188445918 0.188445918
## [385] 0.188445918 0.188445918 0.196368910 0.196368910 0.196368910 0.196368910
## [391] 0.204695686 0.204695686 0.204695686 0.213256293 0.213256293 0.213256293
## [397] 0.213256293 0.213256293 0.221921101 0.221921101 0.230668576 0.239491574
## [403] 0.239491574 0.248425299 0.248425299 0.248425299 0.248425299 0.248425299
## [409] 0.257512175 0.266684841 0.266684841 0.266684841 0.266684841 0.266684841
## [415] 0.266684841 0.275887140 0.275887140 0.285301563 0.285301563 0.294762418
## [421] 0.304319540 0.304319540 0.304319540 0.304319540 0.304319540 0.304319540
## [427] 0.304319540 0.304319540 0.304319540 0.314131316 0.314131316 0.324057279
## [433] 0.334041522 0.334041522 0.344145476 0.344145476 0.344145476 0.344145476
## [439] 0.344145476 0.344145476 0.344145476 0.354489364 0.354489364 0.364955963
## [445] 0.364955963 0.375468287 0.375468287 0.375468287 0.375468287 0.375468287
## [451] 0.386067530 0.396810622 0.396810622 0.396810622 0.418872475 0.418872475
## [457] 0.430217353 0.430217353 0.430217353 0.441771494 0.453445406 0.453445406
## [463] 0.465319253 0.465319253 0.477384731 0.489634257 0.502028814 0.514632492
## [469] 0.527340450 0.540164985 0.540164985 0.540164985 0.540164985 0.553304597
## [475] 0.553304597 0.566657491 0.580311979 0.580311979 0.580311979 0.594305960
## [481] 0.609012152 0.623904225 0.639143157 0.639143157 0.639143157 0.654807088
## [487] 0.670718435 0.670718435 0.687301167 0.704071054 0.704071054 0.704071054
## [493] 0.721157987 0.721157987 0.738561528 0.756346779 0.756346779 0.756346779
## [499] 0.774488197 0.774488197 0.774488197 0.774488197 0.774488197 0.793145557
## [505] 0.793145557 0.812226345 0.812226345 0.812226345 0.831644405 0.851361126
## [511] 0.871448014 0.871448014 0.871448014 0.891810367 0.891810367 0.913268483
## [517] 0.913268483 0.913268483 0.913268483 0.913268483 0.936798339 0.936798339
## [523] 0.936798339 0.936798339 0.936798339 0.961841394 0.987712859 0.987712859
## [529] 1.013960311 1.040953046 1.040953046 1.040953046 1.069941379 1.099543534
## [535] 1.099543534 1.099543534 1.099543534 1.099543534 1.129531930 1.129531930
## [541] 1.160547040 1.192692166 1.228912605 1.228912605 1.228912605 1.228912605
## [547] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [553] 1.228912605 1.270816176 1.314046197 1.314046197 1.359029219 1.404934747
## [559] 1.404934747 1.404934747 1.404934747 1.404934747 1.458863406 1.458863406
## [565] 1.458863406 1.458863406 1.520684154 1.520684154 1.590254427 1.590254427
## [571] 1.667260347 1.667260347 1.747291122 1.747291122 1.747291122 1.747291122
## [577] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [583] 1.747291122 1.747291122 1.747291122 1.926625374 1.926625374 1.926625374
## [589] 1.926625374 1.926625374 2.281783691 2.281783691 0.009633169 0.009633169
## [595] 0.014551755 0.014551755 0.014551755 0.019546218 0.019546218 0.024599344
## [601] 0.029678509 0.034864290 0.040107303 0.045413686 0.045413686 0.050750288
## [607] 0.050750288 0.061572580 0.067117647 0.072742518 0.078415489 0.078415489
## [613] 0.084121378 0.089991681 0.095877420 0.101853535 0.107873484 0.107873484
## [619] 0.107873484 0.113924256 0.119989528 0.119989528 0.126231259 0.126231259
## [625] 0.132564626 0.132564626 0.138970374 0.138970374 0.138970374 0.138970374
## [631] 0.158792804 0.158792804 0.165618521 0.165618521 0.172484282 0.172484282
## [637] 0.186450350 0.193667471 0.200947081 0.208250442 0.223100168 0.230731258
## [643] 0.230731258 0.230731258 0.238451311 0.246241191 0.246241191 0.246241191
## [649] 0.246241191 0.254091780 0.254091780 0.261998834 0.269933185 0.269933185
## [655] 0.295740183 0.295740183 0.304607658 0.304607658 0.304607658 0.313541365
## [661] 0.313541365 0.322590148 0.331761267 0.341155661 0.341155661 0.350635123
## [667] 0.350635123 0.360361864 0.360361864 0.370472540 0.380732678 0.391154668
## [673] 0.391154668 0.401655837 0.401655837 0.401655837 0.412274302 0.412274302
## [679] 0.422970130 0.422970130 0.422970130 0.433912349 0.433912349 0.433912349
## [685] 0.433912349 0.445124076 0.445124076 0.445124076 0.456660219 0.456660219
## [691] 0.456660219 0.468424656 0.480317909 0.480317909 0.480317909 0.492387762
## [697] 0.492387762 0.492387762 0.504593246 0.504593246 0.516909136 0.516909136
## [703] 0.529347066 0.529347066 0.529347066 0.529347066 0.541876644 0.554515396
## [709] 0.567272570 0.580175913 0.593186674 0.593186674 0.606344153 0.606344153
## [715] 0.619675763 0.619675763 0.619675763 0.633344103 0.647161102 0.647161102
## [721] 0.647161102 0.647161102 0.661946165 0.677209373 0.677209373 0.677209373
## [727] 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700
## [733] 0.708979609 0.708979609 0.725379668 0.725379668 0.742673869 0.742673869
## [739] 0.760665334 0.760665334 0.760665334 0.760665334 0.779916835 0.799534319
## [745] 0.819325249 0.839297067 0.839297067 0.839297067 0.859999536 0.859999536
## [751] 0.859999536 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673
## [757] 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673 0.905686340
## [763] 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340
## [769] 0.930593182 0.930593182 0.956024835 0.956024835 0.956024835 0.956024835
## [775] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [781] 0.982860461 0.982860461 0.982860461 1.011138349 1.011138349 1.011138349
## [787] 1.040995345 1.072147856 1.072147856 1.104113664 1.104113664 1.104113664
## [793] 1.136775819 1.171345121 1.171345121 1.207040278 1.207040278 1.207040278
## [799] 1.207040278 1.207040278 1.246907275 1.246907275 1.287930715 1.330411123
## [805] 1.330411123 1.330411123 1.330411123 1.374604413 1.420519943 1.420519943
## [811] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943
## [817] 1.420519943 1.420519943 1.470281401 1.521709021 1.521709021 1.521709021
## [823] 1.521709021 1.574890670 1.574890670 1.574890670 1.574890670 1.574890670
## [829] 1.574890670 1.574890670 1.629630555 1.685791093 1.685791093 1.685791093
## [835] 1.685791093 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [841] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [847] 1.832400858 1.832400858 1.917615102 1.917615102 1.917615102 1.917615102
## [853] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [859] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [865] 1.917615102 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [871] 2.029441753 2.029441753 2.029441753 2.029441753 2.194031226 2.194031226
## [877] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226
## [883] 2.530841239 2.925483659 2.925483659 2.925483659 2.925483659 2.925483659
res4<-residuals(Cox.fit4,type="deviance")
hist(res4)
<<<<<<< HEAD

boxplot(res4~Cdata$Death_Censoring)

=======

boxplot(res4~Cdata$Death_Censoring)

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
# Looking at Protrombin
min(Cdata[,8])# Lowest 12
## [1] 12
max(Cdata[,8])# Highest 135
## [1] 135
Prothrombin.groups=cut(Cdata$Prothrombin,breaks=c(-Inf,quantile(Cdata$Prothrombin ,prob=c(0.33,0.66)),Inf)) #devides the Age into groups
Cox.fit5 <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ strata(Prothrombin.groups), data = Cdata)
summary(Cox.fit5)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + strata(Prothrombin.groups), data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                        coef exp(coef) se(coef)     z Pr(>|z|)    
## factor(Treatment)1 0.100662  1.105902 0.118140 0.852  0.39418    
## factor(Sex)1       0.475374  1.608615 0.126927 3.745  0.00018 ***
## factor(Ascites)1   0.460651  1.585105 0.178184 2.585  0.00973 ** 
## factor(Ascites)2   0.952334  2.591752 0.184248 5.169 2.36e-07 ***
## Age                0.048430  1.049622 0.006825 7.096 1.28e-12 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                    exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1     1.106     0.9042    0.8773     1.394
## factor(Sex)1           1.609     0.6217    1.2543     2.063
## factor(Ascites)1       1.585     0.6309    1.1179     2.248
## factor(Ascites)2       2.592     0.3858    1.8062     3.719
## Age                    1.050     0.9527    1.0357     1.064
## 
## Concordance= 0.672  (se = 0.018 )
## Likelihood ratio test= 91.68  on 5 df,   p=<2e-16
## Wald test            = 90.79  on 5 df,   p=<2e-16
## Score (logrank) test = 93.72  on 5 df,   p=<2e-16
source("andersen.R")
test5<-survfit(Cox.fit5)
test5
## Call: survfit(formula = Cox.fit5)
## 
##             n events median 0.95LCL 0.95UCL
## (-Inf,58] 171    127   1169     883    1582
## (58,80]   165     96   1698    1203    2195
## (80, Inf] 152     69   2364    1976    3132
andersen.plot(test5) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD

=======

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##   [1] 0.003815209 0.007683007 0.011584233 0.015524640 0.019505384 0.023538256
##   [7] 0.023538256 0.027640426 0.035938955 0.035938955 0.040147080 0.044419856
##  [13] 0.044419856 0.044419856 0.044419856 0.044419856 0.048804220 0.062254913
##  [19] 0.066857759 0.076224418 0.076224418 0.076224418 0.076224418 0.080982967
##  [25] 0.085819661 0.090716044 0.095649157 0.095649157 0.100605383 0.100605383
##  [31] 0.105595367 0.105595367 0.110675694 0.115868619 0.115868619 0.115868619
##  [37] 0.121148141 0.121148141 0.121148141 0.121148141 0.121148141 0.131893456
##  [43] 0.131893456 0.137421743 0.142988787 0.142988787 0.148646614 0.148646614
##  [49] 0.160228415 0.160228415 0.166295953 0.178610032 0.178610032 0.191138887
##  [55] 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449
##  [61] 0.204138206 0.210707735 0.217318317 0.224215818 0.238147577 0.238147577
##  [67] 0.245219522 0.252306963 0.259472322 0.259472322 0.274128676 0.274128676
##  [73] 0.274128676 0.281546394 0.281546394 0.289018618 0.289018618 0.296531974
##  [79] 0.304353874 0.312209872 0.312209872 0.320138856 0.320138856 0.320138856
##  [85] 0.328159980 0.336243939 0.336243939 0.344381675 0.344381675 0.352598926
##  [91] 0.352598926 0.360910937 0.369392228 0.378042605 0.378042605 0.378042605
##  [97] 0.378042605 0.386845320 0.386845320 0.395828950 0.395828950 0.395828950
## [103] 0.395828950 0.395828950 0.404948309 0.414138672 0.414138672 0.414138672
## [109] 0.414138672 0.414138672 0.423426979 0.423426979 0.423426979 0.432806013
## [115] 0.432806013 0.442290826 0.442290826 0.451875011 0.451875011 0.461510440
## [121] 0.461510440 0.471258993 0.471258993 0.471258993 0.481150204 0.481150204
## [127] 0.481150204 0.481150204 0.491080246 0.491080246 0.491080246 0.491080246
## [133] 0.491080246 0.501355641 0.511878399 0.511878399 0.511878399 0.522589991
## [139] 0.522589991 0.533493357 0.533493357 0.544539943 0.544539943 0.544539943
## [145] 0.544539943 0.544539943 0.544539943 0.544539943 0.555981864 0.567486954
## [151] 0.567486954 0.579419841 0.579419841 0.579419841 0.579419841 0.579419841
## [157] 0.579419841 0.592453938 0.605747671 0.605747671 0.619209067 0.633357923
## [163] 0.633357923 0.633357923 0.633357923 0.633357923 0.633357923 0.648497735
## [169] 0.663809559 0.663809559 0.663809559 0.663809559 0.663809559 0.679715552
## [175] 0.695849750 0.712319855 0.729088039 0.729088039 0.729088039 0.729088039
## [181] 0.729088039 0.746646274 0.764560852 0.764560852 0.764560852 0.783263260
## [187] 0.802337188 0.802337188 0.822432633 0.842613180 0.842613180 0.863519292
## [193] 0.863519292 0.863519292 0.884708578 0.884708578 0.884708578 0.884708578
## [199] 0.906698043 0.906698043 0.928838016 0.928838016 0.928838016 0.928838016
## [205] 0.928838016 0.951641911 0.951641911 0.951641911 0.951641911 0.951641911
## [211] 0.951641911 0.951641911 0.951641911 0.975097788 0.999170914 0.999170914
## [217] 1.023722913 1.048726373 1.074879725 1.074879725 1.102723976 1.102723976
## [223] 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976
## [229] 1.132454776 1.132454776 1.132454776 1.132454776 1.163457833 1.196856365
## [235] 1.196856365 1.196856365 1.231339476 1.231339476 1.271278329 1.271278329
## [241] 1.271278329 1.315576550 1.315576550 1.360742931 1.360742931 1.360742931
## [247] 1.360742931 1.408744019 1.408744019 1.408744019 1.408744019 1.408744019
## [253] 1.408744019 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [259] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [265] 1.461889660 1.461889660 1.524568595 1.524568595 1.594403830 1.594403830
## [271] 1.668862962 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [277] 1.751026186 1.751026186 1.842255469 1.842255469 1.842255469 1.842255469
## [283] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469
## [289] 1.842255469 1.842255469 1.842255469 1.842255469 2.124490838 2.124490838
## [295] 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161
## [301] 2.496124161 2.496124161 2.496124161 0.000000000 0.000000000 0.000000000
## [307] 0.000000000 0.000000000 0.000000000 0.000000000 0.005214536 0.005214536
## [313] 0.005214536 0.005214536 0.005214536 0.010551942 0.010551942 0.015947912
## [319] 0.015947912 0.015947912 0.015947912 0.015947912 0.015947912 0.021407242
## [325] 0.021407242 0.021407242 0.021407242 0.021407242 0.021407242 0.021407242
## [331] 0.021407242 0.021407242 0.026971836 0.026971836 0.026971836 0.026971836
## [337] 0.026971836 0.032611990 0.032611990 0.032611990 0.038280941 0.038280941
## [343] 0.038280941 0.038280941 0.044119420 0.050049846 0.050049846 0.050049846
## [349] 0.050049846 0.050049846 0.050049846 0.050049846 0.056090678 0.056090678
## [355] 0.056090678 0.062210413 0.062210413 0.062210413 0.062210413 0.068421665
## [361] 0.074704142 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806
## [367] 0.081095806 0.094314691 0.101097016 0.101097016 0.101097016 0.101097016
## [373] 0.101097016 0.101097016 0.107965714 0.115044736 0.115044736 0.122190020
## [379] 0.122190020 0.122190020 0.129556960 0.129556960 0.129556960 0.136977660
## [385] 0.136977660 0.144630677 0.152345180 0.152345180 0.152345180 0.160176278
## [391] 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278
## [397] 0.160176278 0.168197928 0.168197928 0.168197928 0.168197928 0.168197928
## [403] 0.168197928 0.168197928 0.168197928 0.176614537 0.176614537 0.176614537
## [409] 0.176614537 0.185176955 0.193829682 0.193829682 0.202671703 0.211594856
## [415] 0.220626646 0.220626646 0.220626646 0.229793501 0.229793501 0.239077829
## [421] 0.248444048 0.257932939 0.257932939 0.267521287 0.267521287 0.277197004
## [427] 0.286917451 0.286917451 0.286917451 0.296920891 0.307216111 0.307216111
## [433] 0.317632823 0.317632823 0.317632823 0.317632823 0.317632823 0.317632823
## [439] 0.328847214 0.340199006 0.340199006 0.351912350 0.351912350 0.351912350
## [445] 0.351912350 0.351912350 0.363932510 0.363932510 0.376103680 0.388411924
## [451] 0.388411924 0.388411924 0.388411924 0.388411924 0.388411924 0.400827235
## [457] 0.413388661 0.426093822 0.426093822 0.438880811 0.438880811 0.438880811
## [463] 0.451909153 0.451909153 0.465106583 0.478494715 0.492428623 0.492428623
## [469] 0.492428623 0.492428623 0.492428623 0.492428623 0.492428623 0.506710060
## [475] 0.521152134 0.521152134 0.521152134 0.535918656 0.535918656 0.535918656
## [481] 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093
## [487] 0.566225029 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [493] 0.581917832 0.581917832 0.581917832 0.581917832 0.598575995 0.598575995
## [499] 0.615540328 0.633502999 0.651685645 0.651685645 0.670216074 0.670216074
## [505] 0.689056961 0.708270080 0.728416190 0.748931218 0.748931218 0.769827448
## [511] 0.791444671 0.813230253 0.813230253 0.835981763 0.859337551 0.859337551
## [517] 0.859337551 0.859337551 0.883790398 0.883790398 0.883790398 0.883790398
## [523] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [529] 0.910092808 0.937699114 0.966355040 0.966355040 0.995929592 1.026020500
## [535] 1.057173383 1.057173383 1.057173383 1.057173383 1.091579070 1.091579070
## [541] 1.127017733 1.127017733 1.127017733 1.163673474 1.163673474 1.163673474
## [547] 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474
## [553] 1.163673474 1.205390594 1.205390594 1.205390594 1.205390594 1.248789143
## [559] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [565] 1.248789143 1.302049650 1.302049650 1.302049650 1.302049650 1.361380792
## [571] 1.361380792 1.361380792 1.423171756 1.423171756 1.423171756 1.423171756
## [577] 1.423171756 1.492959677 1.565134027 1.641566599 1.641566599 1.641566599
## [583] 1.641566599 1.641566599 1.641566599 1.760594965 1.760594965 1.760594965
## [589] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [595] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.995646102
## [601] 1.995646102 1.995646102 1.995646102 1.995646102 1.995646102 1.995646102
## [607] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [613] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [619] 0.007077094 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278
## [625] 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278
## [631] 0.021477375 0.021477375 0.021477375 0.021477375 0.021477375 0.021477375
## [637] 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193
## [643] 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193 0.043683254
## [649] 0.043683254 0.043683254 0.043683254 0.043683254 0.043683254 0.051238198
## [655] 0.051238198 0.058848646 0.058848646 0.066482062 0.066482062 0.066482062
## [661] 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062
## [667] 0.066482062 0.066482062 0.082161595 0.082161595 0.082161595 0.082161595
## [673] 0.082161595 0.090312041 0.090312041 0.107058768 0.107058768 0.107058768
## [679] 0.107058768 0.107058768 0.115642909 0.115642909 0.115642909 0.124279460
## [685] 0.124279460 0.124279460 0.124279460 0.124279460 0.124279460 0.132971217
## [691] 0.132971217 0.132971217 0.132971217 0.141833373 0.141833373 0.141833373
## [697] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [703] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [709] 0.141833373 0.141833373 0.141833373 0.151197061 0.151197061 0.151197061
## [715] 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061
## [721] 0.160686642 0.160686642 0.170396775 0.170396775 0.180213133 0.190108901
## [727] 0.190108901 0.210270329 0.210270329 0.210270329 0.210270329 0.210270329
## [733] 0.220576216 0.220576216 0.230965174 0.230965174 0.230965174 0.230965174
## [739] 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174
## [745] 0.241881805 0.241881805 0.241881805 0.241881805 0.241881805 0.241881805
## [751] 0.241881805 0.241881805 0.252853192 0.252853192 0.252853192 0.252853192
## [757] 0.252853192 0.252853192 0.252853192 0.252853192 0.252853192 0.264161947
## [763] 0.275576224 0.310613892 0.310613892 0.310613892 0.310613892 0.322632099
## [769] 0.322632099 0.322632099 0.322632099 0.322632099 0.322632099 0.322632099
## [775] 0.322632099 0.322632099 0.334704748 0.334704748 0.334704748 0.334704748
## [781] 0.334704748 0.334704748 0.347024049 0.372137832 0.372137832 0.372137832
## [787] 0.372137832 0.385220199 0.385220199 0.385220199 0.398660095 0.398660095
## [793] 0.398660095 0.398660095 0.412484044 0.426505935 0.426505935 0.426505935
## [799] 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935
## [805] 0.426505935 0.456128884 0.471685153 0.471685153 0.471685153 0.471685153
## [811] 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692
## [817] 0.487657692 0.504636398 0.504636398 0.504636398 0.540398037 0.540398037
## [823] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442
## [829] 0.597092442 0.597092442 0.597092442 0.597092442 0.617805328 0.617805328
## [835] 0.617805328 0.617805328 0.640376113 0.640376113 0.664790137 0.715491227
## [841] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [847] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [853] 0.715491227 0.715491227 0.715491227 0.715491227 0.746629486 0.746629486
## [859] 0.746629486 0.746629486 0.746629486 0.746629486 0.780446351 0.780446351
## [865] 0.851972559 0.851972559 0.934636454 0.934636454 0.934636454 0.934636454
## [871] 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454
## [877] 0.934636454 0.934636454 0.934636454 0.986212450 0.986212450 1.103396433
## [883] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [889] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [895] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.238566987
## [901] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
## [907] 1.390361656 1.390361656 1.390361656
##   [1] 0.003815209 0.007683007 0.011584233 0.015524640 0.019505384 0.023538256
##   [7] 0.023538256 0.023538256 0.027640426 0.035938955 0.035938955 0.040147080
##  [13] 0.044419856 0.044419856 0.044419856 0.044419856 0.044419856 0.048804220
##  [19] 0.062254913 0.066857759 0.076224418 0.076224418 0.080982967 0.080982967
##  [25] 0.085819661 0.090716044 0.095649157 0.100605383 0.100605383 0.100605383
##  [31] 0.105595367 0.110675694 0.115868619 0.115868619 0.121148141 0.121148141
##  [37] 0.131893456 0.137421743 0.142988787 0.142988787 0.148646614 0.148646614
##  [43] 0.160228415 0.160228415 0.166295953 0.178610032 0.191138887 0.197592449
##  [49] 0.197592449 0.197592449 0.197592449 0.204138206 0.210707735 0.210707735
##  [55] 0.210707735 0.210707735 0.217318317 0.224215818 0.238147577 0.245219522
##  [61] 0.252306963 0.252306963 0.259472322 0.259472322 0.259472322 0.274128676
##  [67] 0.274128676 0.281546394 0.281546394 0.289018618 0.289018618 0.296531974
##  [73] 0.304353874 0.312209872 0.320138856 0.320138856 0.320138856 0.328159980
##  [79] 0.336243939 0.344381675 0.352598926 0.352598926 0.360910937 0.369392228
##  [85] 0.378042605 0.386845320 0.395828950 0.395828950 0.404948309 0.404948309
##  [91] 0.404948309 0.414138672 0.414138672 0.414138672 0.423426979 0.432806013
##  [97] 0.432806013 0.442290826 0.442290826 0.442290826 0.442290826 0.442290826
## [103] 0.451875011 0.451875011 0.451875011 0.461510440 0.461510440 0.461510440
## [109] 0.461510440 0.471258993 0.481150204 0.481150204 0.481150204 0.491080246
## [115] 0.491080246 0.491080246 0.501355641 0.511878399 0.511878399 0.522589991
## [121] 0.522589991 0.533493357 0.544539943 0.544539943 0.544539943 0.544539943
## [127] 0.555981864 0.567486954 0.567486954 0.579419841 0.579419841 0.579419841
## [133] 0.579419841 0.579419841 0.579419841 0.579419841 0.592453938 0.605747671
## [139] 0.619209067 0.619209067 0.619209067 0.633357923 0.633357923 0.633357923
## [145] 0.648497735 0.663809559 0.663809559 0.663809559 0.679715552 0.679715552
## [151] 0.695849750 0.712319855 0.712319855 0.712319855 0.729088039 0.729088039
## [157] 0.729088039 0.729088039 0.729088039 0.729088039 0.729088039 0.729088039
## [163] 0.746646274 0.764560852 0.764560852 0.764560852 0.783263260 0.802337188
## [169] 0.802337188 0.822432633 0.822432633 0.822432633 0.842613180 0.842613180
## [175] 0.863519292 0.863519292 0.884708578 0.906698043 0.906698043 0.906698043
## [181] 0.906698043 0.906698043 0.928838016 0.928838016 0.951641911 0.951641911
## [187] 0.951641911 0.951641911 0.975097788 0.975097788 0.975097788 0.975097788
## [193] 0.999170914 0.999170914 0.999170914 0.999170914 1.023722913 1.048726373
## [199] 1.074879725 1.074879725 1.102723976 1.102723976 1.102723976 1.102723976
## [205] 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976
## [211] 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776
## [217] 1.132454776 1.163457833 1.163457833 1.163457833 1.196856365 1.196856365
## [223] 1.196856365 1.231339476 1.271278329 1.271278329 1.315576550 1.360742931
## [229] 1.360742931 1.360742931 1.360742931 1.360742931 1.360742931 1.360742931
## [235] 1.408744019 1.408744019 1.408744019 1.408744019 1.408744019 1.408744019
## [241] 1.408744019 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [247] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [253] 1.461889660 1.524568595 1.524568595 1.594403830 1.594403830 1.668862962
## [259] 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [265] 1.751026186 1.751026186 1.751026186 1.751026186 1.842255469 1.842255469
## [271] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469
## [277] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 2.124490838
## [283] 2.124490838 2.124490838 2.124490838 2.496124161 2.496124161 2.496124161
## [289] 2.496124161 2.496124161 2.496124161 2.496124161 0.000000000 0.000000000
## [295] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [301] 0.005214536 0.005214536 0.005214536 0.005214536 0.005214536 0.010551942
## [307] 0.010551942 0.010551942 0.015947912 0.015947912 0.015947912 0.015947912
## [313] 0.015947912 0.021407242 0.021407242 0.021407242 0.021407242 0.021407242
## [319] 0.021407242 0.021407242 0.026971836 0.026971836 0.026971836 0.026971836
## [325] 0.026971836 0.032611990 0.032611990 0.038280941 0.044119420 0.050049846
## [331] 0.050049846 0.050049846 0.050049846 0.050049846 0.050049846 0.050049846
## [337] 0.056090678 0.056090678 0.062210413 0.062210413 0.081095806 0.081095806
## [343] 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806
## [349] 0.081095806 0.081095806 0.094314691 0.101097016 0.101097016 0.101097016
## [355] 0.101097016 0.101097016 0.101097016 0.101097016 0.115044736 0.115044736
## [361] 0.122190020 0.122190020 0.122190020 0.129556960 0.129556960 0.129556960
## [367] 0.136977660 0.144630677 0.152345180 0.152345180 0.152345180 0.160176278
## [373] 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278 0.168197928
## [379] 0.168197928 0.168197928 0.176614537 0.176614537 0.176614537 0.176614537
## [385] 0.176614537 0.202671703 0.211594856 0.220626646 0.220626646 0.229793501
## [391] 0.229793501 0.229793501 0.229793501 0.239077829 0.248444048 0.248444048
## [397] 0.257932939 0.257932939 0.267521287 0.267521287 0.267521287 0.267521287
## [403] 0.286917451 0.286917451 0.307216111 0.307216111 0.317632823 0.317632823
## [409] 0.317632823 0.317632823 0.317632823 0.340199006 0.340199006 0.351912350
## [415] 0.351912350 0.363932510 0.376103680 0.388411924 0.388411924 0.388411924
## [421] 0.388411924 0.388411924 0.426093822 0.426093822 0.426093822 0.438880811
## [427] 0.438880811 0.438880811 0.438880811 0.438880811 0.451909153 0.451909153
## [433] 0.451909153 0.465106583 0.492428623 0.492428623 0.492428623 0.492428623
## [439] 0.492428623 0.492428623 0.521152134 0.521152134 0.535918656 0.535918656
## [445] 0.535918656 0.535918656 0.535918656 0.535918656 0.535918656 0.550949093
## [451] 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093
## [457] 0.566225029 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [463] 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [469] 0.598575995 0.651685645 0.651685645 0.651685645 0.670216074 0.670216074
## [475] 0.670216074 0.728416190 0.748931218 0.813230253 0.813230253 0.835981763
## [481] 0.859337551 0.859337551 0.859337551 0.859337551 0.859337551 0.883790398
## [487] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [493] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [499] 0.910092808 0.910092808 0.966355040 0.966355040 0.966355040 0.995929592
## [505] 0.995929592 0.995929592 1.057173383 1.057173383 1.057173383 1.057173383
## [511] 1.057173383 1.057173383 1.057173383 1.057173383 1.091579070 1.091579070
## [517] 1.127017733 1.127017733 1.163673474 1.163673474 1.163673474 1.163673474
## [523] 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474
## [529] 1.163673474 1.163673474 1.163673474 1.205390594 1.205390594 1.205390594
## [535] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [541] 1.248789143 1.248789143 1.248789143 1.248789143 1.302049650 1.361380792
## [547] 1.361380792 1.361380792 1.423171756 1.423171756 1.423171756 1.423171756
## [553] 1.423171756 1.423171756 1.492959677 1.492959677 1.492959677 1.641566599
## [559] 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599
## [565] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [571] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [577] 1.760594965 1.760594965 1.760594965 1.760594965 1.995646102 1.995646102
## [583] 1.995646102 1.995646102 0.000000000 0.000000000 0.000000000 0.000000000
## [589] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [595] 0.000000000 0.000000000 0.000000000 0.007077094 0.007077094 0.014216278
## [601] 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278
## [607] 0.014216278 0.021477375 0.021477375 0.021477375 0.021477375 0.021477375
## [613] 0.028820277 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193
## [619] 0.036190193 0.036190193 0.043683254 0.043683254 0.043683254 0.043683254
## [625] 0.043683254 0.051238198 0.051238198 0.058848646 0.058848646 0.066482062
## [631] 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062
## [637] 0.066482062 0.066482062 0.074283328 0.082161595 0.082161595 0.082161595
## [643] 0.082161595 0.082161595 0.090312041 0.090312041 0.090312041 0.098605212
## [649] 0.107058768 0.107058768 0.107058768 0.107058768 0.115642909 0.115642909
## [655] 0.124279460 0.124279460 0.124279460 0.124279460 0.124279460 0.132971217
## [661] 0.132971217 0.132971217 0.132971217 0.141833373 0.141833373 0.141833373
## [667] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [673] 0.141833373 0.141833373 0.151197061 0.151197061 0.151197061 0.151197061
## [679] 0.151197061 0.151197061 0.160686642 0.160686642 0.160686642 0.160686642
## [685] 0.170396775 0.170396775 0.170396775 0.180213133 0.190108901 0.190108901
## [691] 0.200148605 0.210270329 0.210270329 0.210270329 0.210270329 0.220576216
## [697] 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174
## [703] 0.230965174 0.230965174 0.241881805 0.241881805 0.241881805 0.241881805
## [709] 0.252853192 0.252853192 0.252853192 0.252853192 0.252853192 0.252853192
## [715] 0.264161947 0.264161947 0.275576224 0.287173209 0.298818473 0.310613892
## [721] 0.310613892 0.310613892 0.310613892 0.310613892 0.322632099 0.322632099
## [727] 0.322632099 0.322632099 0.322632099 0.322632099 0.322632099 0.334704748
## [733] 0.334704748 0.334704748 0.334704748 0.334704748 0.334704748 0.347024049
## [739] 0.347024049 0.359487200 0.372137832 0.372137832 0.372137832 0.372137832
## [745] 0.372137832 0.385220199 0.385220199 0.385220199 0.398660095 0.398660095
## [751] 0.398660095 0.398660095 0.412484044 0.412484044 0.412484044 0.426505935
## [757] 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935
## [763] 0.441039991 0.456128884 0.456128884 0.471685153 0.471685153 0.487657692
## [769] 0.487657692 0.487657692 0.487657692 0.504636398 0.504636398 0.504636398
## [775] 0.522274168 0.540398037 0.540398037 0.558903304 0.577733813 0.597092442
## [781] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442
## [787] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.617805328
## [793] 0.617805328 0.617805328 0.617805328 0.617805328 0.617805328 0.640376113
## [799] 0.640376113 0.640376113 0.664790137 0.664790137 0.689971460 0.715491227
## [805] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [811] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [817] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [823] 0.746629486 0.746629486 0.746629486 0.746629486 0.780446351 0.815837029
## [829] 0.851972559 0.851972559 0.851972559 0.851972559 0.892199662 0.934636454
## [835] 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454
## [841] 0.934636454 0.934636454 0.934636454 0.934636454 0.986212450 0.986212450
## [847] 1.043419749 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [853] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [859] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [865] 1.103396433 1.103396433 1.238566987 1.238566987 1.390361656 1.390361656
## [871] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
##   [1] 0.023538256 0.027640426 0.035938955 0.035938955 0.044419856 0.044419856
##   [7] 0.044419856 0.044419856 0.044419856 0.076224418 0.076224418 0.080982967
##  [13] 0.095649157 0.100605383 0.100605383 0.100605383 0.105595367 0.110675694
##  [19] 0.115868619 0.121148141 0.121148141 0.121148141 0.121148141 0.131893456
##  [25] 0.131893456 0.148646614 0.160228415 0.160228415 0.178610032 0.178610032
##  [31] 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449
##  [37] 0.210707735 0.210707735 0.210707735 0.217318317 0.224215818 0.238147577
##  [43] 0.238147577 0.252306963 0.252306963 0.259472322 0.259472322 0.259472322
##  [49] 0.274128676 0.274128676 0.274128676 0.281546394 0.289018618 0.289018618
##  [55] 0.296531974 0.312209872 0.320138856 0.320138856 0.320138856 0.336243939
##  [61] 0.344381675 0.344381675 0.352598926 0.352598926 0.378042605 0.378042605
##  [67] 0.378042605 0.386845320 0.395828950 0.395828950 0.395828950 0.395828950
##  [73] 0.395828950 0.404948309 0.404948309 0.414138672 0.414138672 0.414138672
##  [79] 0.414138672 0.414138672 0.414138672 0.423426979 0.423426979 0.423426979
##  [85] 0.432806013 0.432806013 0.442290826 0.442290826 0.442290826 0.442290826
##  [91] 0.442290826 0.451875011 0.451875011 0.451875011 0.451875011 0.461510440
##  [97] 0.461510440 0.461510440 0.461510440 0.471258993 0.471258993 0.481150204
## [103] 0.481150204 0.481150204 0.481150204 0.481150204 0.491080246 0.491080246
## [109] 0.491080246 0.491080246 0.501355641 0.511878399 0.511878399 0.511878399
## [115] 0.522589991 0.522589991 0.533493357 0.544539943 0.544539943 0.544539943
## [121] 0.544539943 0.544539943 0.579419841 0.579419841 0.579419841 0.579419841
## [127] 0.579419841 0.579419841 0.579419841 0.579419841 0.579419841 0.605747671
## [133] 0.619209067 0.619209067 0.633357923 0.633357923 0.633357923 0.633357923
## [139] 0.663809559 0.663809559 0.663809559 0.663809559 0.679715552 0.695849750
## [145] 0.712319855 0.712319855 0.729088039 0.729088039 0.729088039 0.729088039
## [151] 0.729088039 0.729088039 0.764560852 0.764560852 0.764560852 0.764560852
## [157] 0.783263260 0.802337188 0.802337188 0.822432633 0.822432633 0.863519292
## [163] 0.863519292 0.863519292 0.884708578 0.884708578 0.884708578 0.906698043
## [169] 0.906698043 0.906698043 0.906698043 0.906698043 0.928838016 0.928838016
## [175] 0.928838016 0.928838016 0.928838016 0.951641911 0.951641911 0.951641911
## [181] 0.951641911 0.951641911 0.951641911 0.951641911 0.951641911 0.951641911
## [187] 0.951641911 0.975097788 0.975097788 0.975097788 0.999170914 0.999170914
## [193] 0.999170914 0.999170914 1.102723976 1.102723976 1.102723976 1.102723976
## [199] 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976
## [205] 1.102723976 1.102723976 1.132454776 1.132454776 1.132454776 1.132454776
## [211] 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776
## [217] 1.163457833 1.163457833 1.196856365 1.196856365 1.196856365 1.196856365
## [223] 1.231339476 1.271278329 1.315576550 1.360742931 1.360742931 1.360742931
## [229] 1.360742931 1.360742931 1.360742931 1.360742931 1.408744019 1.408744019
## [235] 1.408744019 1.408744019 1.408744019 1.461889660 1.461889660 1.461889660
## [241] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [247] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [253] 1.461889660 1.461889660 1.461889660 1.461889660 1.524568595 1.594403830
## [259] 1.594403830 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [265] 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [271] 1.751026186 1.751026186 1.842255469 1.842255469 1.842255469 1.842255469
## [277] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469
## [283] 1.842255469 2.124490838 2.124490838 2.496124161 2.496124161 2.496124161
## [289] 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161
## [295] 2.496124161 0.000000000 0.005214536 0.005214536 0.005214536 0.010551942
## [301] 0.010551942 0.010551942 0.015947912 0.015947912 0.021407242 0.021407242
## [307] 0.021407242 0.021407242 0.026971836 0.026971836 0.026971836 0.026971836
## [313] 0.026971836 0.032611990 0.038280941 0.038280941 0.038280941 0.038280941
## [319] 0.044119420 0.050049846 0.050049846 0.050049846 0.056090678 0.056090678
## [325] 0.062210413 0.062210413 0.068421665 0.074704142 0.081095806 0.081095806
## [331] 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806
## [337] 0.094314691 0.101097016 0.101097016 0.101097016 0.101097016 0.101097016
## [343] 0.101097016 0.107965714 0.115044736 0.115044736 0.122190020 0.122190020
## [349] 0.122190020 0.129556960 0.136977660 0.144630677 0.144630677 0.152345180
## [355] 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278 0.168197928
## [361] 0.168197928 0.168197928 0.168197928 0.168197928 0.168197928 0.168197928
## [367] 0.176614537 0.176614537 0.176614537 0.176614537 0.176614537 0.185176955
## [373] 0.193829682 0.193829682 0.202671703 0.202671703 0.211594856 0.220626646
## [379] 0.220626646 0.220626646 0.229793501 0.229793501 0.229793501 0.229793501
## [385] 0.239077829 0.239077829 0.248444048 0.248444048 0.257932939 0.257932939
## [391] 0.267521287 0.267521287 0.267521287 0.267521287 0.277197004 0.286917451
## [397] 0.286917451 0.286917451 0.296920891 0.307216111 0.307216111 0.317632823
## [403] 0.317632823 0.317632823 0.317632823 0.317632823 0.317632823 0.328847214
## [409] 0.340199006 0.340199006 0.351912350 0.351912350 0.351912350 0.363932510
## [415] 0.376103680 0.376103680 0.388411924 0.400827235 0.413388661 0.426093822
## [421] 0.426093822 0.426093822 0.438880811 0.438880811 0.438880811 0.438880811
## [427] 0.451909153 0.451909153 0.451909153 0.465106583 0.478494715 0.492428623
## [433] 0.492428623 0.492428623 0.506710060 0.521152134 0.521152134 0.521152134
## [439] 0.535918656 0.535918656 0.535918656 0.535918656 0.535918656 0.550949093
## [445] 0.550949093 0.550949093 0.550949093 0.566225029 0.566225029 0.581917832
## [451] 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [457] 0.581917832 0.581917832 0.598575995 0.615540328 0.633502999 0.651685645
## [463] 0.651685645 0.651685645 0.670216074 0.670216074 0.670216074 0.689056961
## [469] 0.708270080 0.728416190 0.728416190 0.748931218 0.769827448 0.791444671
## [475] 0.813230253 0.813230253 0.813230253 0.813230253 0.835981763 0.835981763
## [481] 0.859337551 0.859337551 0.859337551 0.859337551 0.859337551 0.883790398
## [487] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [493] 0.883790398 0.883790398 0.910092808 0.910092808 0.910092808 0.937699114
## [499] 0.966355040 0.966355040 0.966355040 0.966355040 0.995929592 0.995929592
## [505] 0.995929592 0.995929592 1.026020500 1.057173383 1.057173383 1.057173383
## [511] 1.057173383 1.057173383 1.057173383 1.057173383 1.057173383 1.091579070
## [517] 1.091579070 1.127017733 1.163673474 1.163673474 1.163673474 1.163673474
## [523] 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474
## [529] 1.163673474 1.205390594 1.205390594 1.205390594 1.248789143 1.248789143
## [535] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [541] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [547] 1.248789143 1.248789143 1.302049650 1.302049650 1.302049650 1.361380792
## [553] 1.423171756 1.423171756 1.423171756 1.423171756 1.423171756 1.492959677
## [559] 1.492959677 1.492959677 1.492959677 1.565134027 1.641566599 1.641566599
## [565] 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599
## [571] 1.641566599 1.641566599 1.760594965 1.760594965 1.760594965 1.760594965
## [577] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [583] 1.760594965 1.995646102 1.995646102 1.995646102 1.995646102 1.995646102
## [589] 1.995646102 1.995646102 1.995646102 0.000000000 0.000000000 0.000000000
## [595] 0.000000000 0.007077094 0.007077094 0.014216278 0.014216278 0.014216278
## [601] 0.014216278 0.014216278 0.021477375 0.021477375 0.021477375 0.028820277
## [607] 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193
## [613] 0.036190193 0.036190193 0.043683254 0.043683254 0.051238198 0.058848646
## [619] 0.058848646 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062
## [625] 0.066482062 0.066482062 0.066482062 0.066482062 0.074283328 0.082161595
## [631] 0.082161595 0.082161595 0.082161595 0.082161595 0.090312041 0.090312041
## [637] 0.098605212 0.107058768 0.107058768 0.107058768 0.107058768 0.107058768
## [643] 0.115642909 0.115642909 0.124279460 0.124279460 0.124279460 0.124279460
## [649] 0.132971217 0.132971217 0.132971217 0.141833373 0.141833373 0.141833373
## [655] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [661] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.151197061
## [667] 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061
## [673] 0.151197061 0.151197061 0.151197061 0.160686642 0.160686642 0.160686642
## [679] 0.160686642 0.170396775 0.170396775 0.170396775 0.170396775 0.180213133
## [685] 0.180213133 0.190108901 0.190108901 0.200148605 0.210270329 0.210270329
## [691] 0.210270329 0.210270329 0.210270329 0.220576216 0.220576216 0.220576216
## [697] 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174
## [703] 0.230965174 0.230965174 0.230965174 0.241881805 0.241881805 0.241881805
## [709] 0.241881805 0.241881805 0.241881805 0.252853192 0.252853192 0.252853192
## [715] 0.252853192 0.252853192 0.264161947 0.275576224 0.275576224 0.287173209
## [721] 0.298818473 0.310613892 0.310613892 0.310613892 0.322632099 0.322632099
## [727] 0.322632099 0.322632099 0.322632099 0.334704748 0.334704748 0.334704748
## [733] 0.334704748 0.334704748 0.334704748 0.334704748 0.347024049 0.359487200
## [739] 0.372137832 0.372137832 0.372137832 0.372137832 0.385220199 0.385220199
## [745] 0.398660095 0.398660095 0.398660095 0.398660095 0.398660095 0.412484044
## [751] 0.412484044 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935
## [757] 0.426505935 0.426505935 0.441039991 0.456128884 0.456128884 0.456128884
## [763] 0.471685153 0.471685153 0.471685153 0.471685153 0.487657692 0.487657692
## [769] 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692
## [775] 0.487657692 0.504636398 0.504636398 0.504636398 0.504636398 0.522274168
## [781] 0.540398037 0.540398037 0.558903304 0.577733813 0.597092442 0.597092442
## [787] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442
## [793] 0.617805328 0.617805328 0.617805328 0.617805328 0.617805328 0.617805328
## [799] 0.617805328 0.617805328 0.617805328 0.640376113 0.640376113 0.640376113
## [805] 0.640376113 0.640376113 0.664790137 0.689971460 0.715491227 0.715491227
## [811] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [817] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [823] 0.715491227 0.715491227 0.746629486 0.746629486 0.746629486 0.746629486
## [829] 0.746629486 0.746629486 0.780446351 0.780446351 0.780446351 0.815837029
## [835] 0.851972559 0.851972559 0.851972559 0.851972559 0.892199662 0.934636454
## [841] 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454
## [847] 0.934636454 0.934636454 0.934636454 0.934636454 0.986212450 0.986212450
## [853] 0.986212450 0.986212450 1.043419749 1.103396433 1.103396433 1.103396433
## [859] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [865] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [871] 1.103396433 1.103396433 1.103396433 1.103396433 1.238566987 1.390361656
## [877] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
## [883] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
res5<-residuals(Cox.fit5,type="deviance")
hist(res5)
<<<<<<< HEAD

boxplot(res5~Cdata$Death_Censoring)

Opgave 3

=======

boxplot(res5~Cdata$Death_Censoring)

Opgave 3

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##                         coef exp(coef)  se(coef)      z        p
## factor(Treatment)1  0.067070  1.069371  0.117868  0.569   0.5693
## factor(Sex)1        0.525289  1.690948  0.126974  4.137 3.52e-05
## factor(Ascites)1    0.441588  1.555174  0.179236  2.464   0.0138
## factor(Ascites)2    0.974063  2.648684  0.181719  5.360 8.31e-08
## Age                 0.049151  1.050378  0.006764  7.266 3.70e-13
## Prothrombin        -0.012923  0.987161  0.002917 -4.431 9.39e-06
## 
## Likelihood ratio test=129.5  on 6 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                         coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Treatment)1  0.067070  1.069371  0.117868  0.569   0.5693    
## factor(Sex)1        0.525289  1.690948  0.126974  4.137 3.52e-05 ***
## factor(Ascites)1    0.441588  1.555174  0.179236  2.464   0.0138 *  
## factor(Ascites)2    0.974063  2.648684  0.181719  5.360 8.31e-08 ***
## Age                 0.049151  1.050378  0.006764  7.266 3.70e-13 ***
## Prothrombin        -0.012923  0.987161  0.002917 -4.431 9.39e-06 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                    exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1    1.0694     0.9351    0.8488    1.3473
## factor(Sex)1          1.6909     0.5914    1.3184    2.1688
## factor(Ascites)1      1.5552     0.6430    1.0945    2.2098
## factor(Ascites)2      2.6487     0.3775    1.8550    3.7819
## Age                   1.0504     0.9520    1.0365    1.0644
## Prothrombin           0.9872     1.0130    0.9815    0.9928
## 
## Concordance= 0.696  (se = 0.017 )
## Likelihood ratio test= 129.5  on 6 df,   p=<2e-16
## Wald test            = 132.3  on 6 df,   p=<2e-16
## Score (logrank) test = 142.4  on 6 df,   p=<2e-16
<<<<<<< HEAD
temp=cox.zph(Cox.fit)
par(mfrow=c(1,2))
plot(temp)

Opgave 4

=======

Opgave 4

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
# pull out different variables to see the effect, leveled varables are not of intresed since all information have be extracted?
# Looking at the model without Age
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit._age <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Prothrombin, data = Cdata)
res_age=residuals(Cox.fit._age)
scatter.smooth(Cdata$Age, res_age)
<<<<<<< HEAD

=======

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
#Looking at the model without Prothrombin
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit._Prothrombin <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age, data = Cdata)
res_Prothrombin=residuals(Cox.fit._Prothrombin)
scatter.smooth(Cdata$Prothrombin, res_Prothrombin)
<<<<<<< HEAD

Opgave 5

=======

Opgave 5

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
#Interaction Treatment and Sex
Cox.fit_TS <- coxph(Sur1~ factor(Treatment)*factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit_TS
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) * factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##                                      coef exp(coef)  se(coef)      z        p
## factor(Treatment)1               0.220814  1.247091  0.192829  1.145 0.252157
## factor(Sex)1                     0.653601  1.922450  0.181024  3.611 0.000306
## factor(Ascites)1                 0.443838  1.558678  0.179358  2.475 0.013339
## factor(Ascites)2                 0.967216  2.630611  0.182275  5.306 1.12e-07
## Age                              0.049307  1.050543  0.006805  7.245 4.32e-13
## Prothrombin                     -0.012752  0.987329  0.002927 -4.357 1.32e-05
## factor(Treatment)1:factor(Sex)1 -0.246543  0.781497  0.243803 -1.011 0.311901
## 
## Likelihood ratio test=130.6  on 7 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TS)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) * factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                                      coef exp(coef)  se(coef)      z Pr(>|z|)
## factor(Treatment)1               0.220814  1.247091  0.192829  1.145 0.252157
## factor(Sex)1                     0.653601  1.922450  0.181024  3.611 0.000306
## factor(Ascites)1                 0.443838  1.558678  0.179358  2.475 0.013339
## factor(Ascites)2                 0.967216  2.630611  0.182275  5.306 1.12e-07
## Age                              0.049307  1.050543  0.006805  7.245 4.32e-13
## Prothrombin                     -0.012752  0.987329  0.002927 -4.357 1.32e-05
## factor(Treatment)1:factor(Sex)1 -0.246543  0.781497  0.243803 -1.011 0.311901
##                                    
## factor(Treatment)1                 
## factor(Sex)1                    ***
## factor(Ascites)1                *  
## factor(Ascites)2                ***
## Age                             ***
## Prothrombin                     ***
## factor(Treatment)1:factor(Sex)1    
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                                 exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1                 1.2471     0.8019    0.8546     1.820
## factor(Sex)1                       1.9225     0.5202    1.3482     2.741
## factor(Ascites)1                   1.5587     0.6416    1.0967     2.215
## factor(Ascites)2                   2.6306     0.3801    1.8404     3.760
## Age                                1.0505     0.9519    1.0366     1.065
## Prothrombin                        0.9873     1.0128    0.9817     0.993
## factor(Treatment)1:factor(Sex)1    0.7815     1.2796    0.4846     1.260
## 
## Concordance= 0.695  (se = 0.017 )
## Likelihood ratio test= 130.6  on 7 df,   p=<2e-16
## Wald test            = 132.1  on 7 df,   p=<2e-16
## Score (logrank) test = 142.7  on 7 df,   p=<2e-16
#Interaction Treatment and Ascites
Cox.fit_TA <- coxph(Sur1~ factor(Sex)+factor(Treatment)*factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit_TA
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Treatment) * factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##                                          coef exp(coef)  se(coef)      z
## factor(Sex)1                         0.503528  1.654548  0.126989  3.965
## factor(Treatment)1                   0.291419  1.338326  0.139323  2.092
## factor(Ascites)1                     0.786840  2.196445  0.234027  3.362
## factor(Ascites)2                     1.455465  4.286476  0.243010  5.989
## Age                                  0.046595  1.047697  0.006718  6.936
## Prothrombin                         -0.013059  0.987025  0.002968 -4.400
## factor(Treatment)1:factor(Ascites)1 -0.729933  0.481941  0.360030 -2.027
## factor(Treatment)1:factor(Ascites)2 -0.902858  0.405409  0.345164 -2.616
##                                            p
## factor(Sex)1                        7.34e-05
## factor(Treatment)1                  0.036467
## factor(Ascites)1                    0.000773
## factor(Ascites)2                    2.11e-09
## Age                                 4.03e-12
## Prothrombin                         1.08e-05
## factor(Treatment)1:factor(Ascites)1 0.042619
## factor(Treatment)1:factor(Ascites)2 0.008904
## 
## Likelihood ratio test=139.1  on 8 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TA)
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Treatment) * factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                                          coef exp(coef)  se(coef)      z
## factor(Sex)1                         0.503528  1.654548  0.126989  3.965
## factor(Treatment)1                   0.291419  1.338326  0.139323  2.092
## factor(Ascites)1                     0.786840  2.196445  0.234027  3.362
## factor(Ascites)2                     1.455465  4.286476  0.243010  5.989
## Age                                  0.046595  1.047697  0.006718  6.936
## Prothrombin                         -0.013059  0.987025  0.002968 -4.400
## factor(Treatment)1:factor(Ascites)1 -0.729933  0.481941  0.360030 -2.027
## factor(Treatment)1:factor(Ascites)2 -0.902858  0.405409  0.345164 -2.616
##                                     Pr(>|z|)    
## factor(Sex)1                        7.34e-05 ***
## factor(Treatment)1                  0.036467 *  
## factor(Ascites)1                    0.000773 ***
## factor(Ascites)2                    2.11e-09 ***
## Age                                 4.03e-12 ***
## Prothrombin                         1.08e-05 ***
## factor(Treatment)1:factor(Ascites)1 0.042619 *  
## factor(Treatment)1:factor(Ascites)2 0.008904 ** 
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                                     exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1                           1.6545     0.6044    1.2900    2.1221
## factor(Treatment)1                     1.3383     0.7472    1.0185    1.7585
## factor(Ascites)1                       2.1964     0.4553    1.3884    3.4748
## factor(Ascites)2                       4.2865     0.2333    2.6623    6.9016
## Age                                    1.0477     0.9545    1.0340    1.0616
## Prothrombin                            0.9870     1.0131    0.9813    0.9928
## factor(Treatment)1:factor(Ascites)1    0.4819     2.0749    0.2380    0.9760
## factor(Treatment)1:factor(Ascites)2    0.4054     2.4666    0.2061    0.7974
## 
## Concordance= 0.701  (se = 0.017 )
## Likelihood ratio test= 139.1  on 8 df,   p=<2e-16
## Wald test            = 148.8  on 8 df,   p=<2e-16
## Score (logrank) test = 167.7  on 8 df,   p=<2e-16
#standard model
Cox.fit <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##                         coef exp(coef)  se(coef)      z        p
## factor(Treatment)1  0.067070  1.069371  0.117868  0.569   0.5693
## factor(Sex)1        0.525289  1.690948  0.126974  4.137 3.52e-05
## factor(Ascites)1    0.441588  1.555174  0.179236  2.464   0.0138
## factor(Ascites)2    0.974063  2.648684  0.181719  5.360 8.31e-08
## Age                 0.049151  1.050378  0.006764  7.266 3.70e-13
## Prothrombin        -0.012923  0.987161  0.002917 -4.431 9.39e-06
## 
## Likelihood ratio test=129.5  on 6 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                         coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Treatment)1  0.067070  1.069371  0.117868  0.569   0.5693    
## factor(Sex)1        0.525289  1.690948  0.126974  4.137 3.52e-05 ***
## factor(Ascites)1    0.441588  1.555174  0.179236  2.464   0.0138 *  
## factor(Ascites)2    0.974063  2.648684  0.181719  5.360 8.31e-08 ***
## Age                 0.049151  1.050378  0.006764  7.266 3.70e-13 ***
## Prothrombin        -0.012923  0.987161  0.002917 -4.431 9.39e-06 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                    exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1    1.0694     0.9351    0.8488    1.3473
## factor(Sex)1          1.6909     0.5914    1.3184    2.1688
## factor(Ascites)1      1.5552     0.6430    1.0945    2.2098
## factor(Ascites)2      2.6487     0.3775    1.8550    3.7819
## Age                   1.0504     0.9520    1.0365    1.0644
## Prothrombin           0.9872     1.0130    0.9815    0.9928
## 
## Concordance= 0.696  (se = 0.017 )
## Likelihood ratio test= 129.5  on 6 df,   p=<2e-16
## Wald test            = 132.3  on 6 df,   p=<2e-16
## Score (logrank) test = 142.4  on 6 df,   p=<2e-16
<<<<<<< HEAD
# note here that Ascites contains 3 variables, so the p-values can't be interpreted directly so we compair Likelihood ratio tests against the standard model
=======
# note here that Ascites contains 3 variables, so the p-values can't be interpreted directly so we compair Likelihood ratio tests against the standard model
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
#we use a anova co get the chi squared value for the 2 models, the difference in likelihood ratio test is 139.1-129.5=9.4 and the DF is 8-6=2
anova(Cox.fit_TA,Cox.fit) # we get at p-value of 0.008 so it is significant.
## Analysis of Deviance Table
##  Cox model: response is  Sur1
##  Model 1: ~ factor(Sex) + factor(Treatment) * factor(Ascites) + Age + Prothrombin
##  Model 2: ~ factor(Treatment) + factor(Sex) + factor(Ascites) + Age + Prothrombin
##    loglik Chisq Df P(>|Chi|)   
## 1 -1519.4                      
## 2 -1524.2 9.558  2  0.008404 **
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
======= ## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
#Interaction Treatment and Age
Cox.fit_TAA <- coxph(Sur1~ factor(Sex)+factor(Ascites)+ factor(Treatment)*Age+ Prothrombin, data = Cdata)
Cox.fit_TAA
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + factor(Treatment) * 
##     Age + Prothrombin, data = Cdata)
## 
##                             coef exp(coef)  se(coef)      z        p
## factor(Sex)1            0.529397  1.697909  0.126685  4.179 2.93e-05
## factor(Ascites)1        0.462584  1.588172  0.180018  2.570   0.0102
## factor(Ascites)2        0.992521  2.698027  0.182583  5.436 5.45e-08
## factor(Treatment)1     -0.964723  0.381089  0.805313 -1.198   0.2309
## Age                     0.040708  1.041548  0.009283  4.385 1.16e-05
## Prothrombin            -0.013056  0.987029  0.002919 -4.473 7.71e-06
## factor(Treatment)1:Age  0.016710  1.016851  0.012905  1.295   0.1954
## 
## Likelihood ratio test=131.2  on 7 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TAA)
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + factor(Treatment) * 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                             coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Sex)1            0.529397  1.697909  0.126685  4.179 2.93e-05 ***
## factor(Ascites)1        0.462584  1.588172  0.180018  2.570   0.0102 *  
## factor(Ascites)2        0.992521  2.698027  0.182583  5.436 5.45e-08 ***
## factor(Treatment)1     -0.964723  0.381089  0.805313 -1.198   0.2309    
## Age                     0.040708  1.041548  0.009283  4.385 1.16e-05 ***
## Prothrombin            -0.013056  0.987029  0.002919 -4.473 7.71e-06 ***
## factor(Treatment)1:Age  0.016710  1.016851  0.012905  1.295   0.1954    
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                        exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1              1.6979     0.5890   1.32458    2.1765
## factor(Ascites)1          1.5882     0.6297   1.11600    2.2601
## factor(Ascites)2          2.6980     0.3706   1.88639    3.8589
## factor(Treatment)1        0.3811     2.6241   0.07862    1.8472
## Age                       1.0415     0.9601   1.02277    1.0607
## Prothrombin               0.9870     1.0131   0.98140    0.9927
## factor(Treatment)1:Age    1.0169     0.9834   0.99145    1.0429
## 
## Concordance= 0.695  (se = 0.017 )
## Likelihood ratio test= 131.2  on 7 df,   p=<2e-16
## Wald test            = 134.3  on 7 df,   p=<2e-16
## Score (logrank) test = 144.7  on 7 df,   p=<2e-16
#Interaction Treatment and Prothrombin
Cox.fit_TP <- coxph(Sur1~ factor(Sex)+factor(Ascites)+ Age+factor(Treatment)*Prothrombin, data = Cdata)
Cox.fit_TP
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + Age + 
##     factor(Treatment) * Prothrombin, data = Cdata)
## 
##                                     coef exp(coef)  se(coef)      z        p
## factor(Sex)1                    0.518523  1.679545  0.127462  4.068 4.74e-05
## factor(Ascites)1                0.441740  1.555411  0.179070  2.467 0.013631
## factor(Ascites)2                0.978996  2.661782  0.181549  5.392 6.95e-08
## Age                             0.048780  1.049989  0.006773  7.202 5.93e-13
## factor(Treatment)1             -0.157295  0.854452  0.377810 -0.416 0.677166
## Prothrombin                    -0.014688  0.985419  0.004080 -3.600 0.000318
## factor(Treatment)1:Prothrombin  0.003450  1.003456  0.005522  0.625 0.532106
## 
## Likelihood ratio test=129.9  on 7 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TP)
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + Age + 
##     factor(Treatment) * Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                                     coef exp(coef)  se(coef)      z Pr(>|z|)
## factor(Sex)1                    0.518523  1.679545  0.127462  4.068 4.74e-05
## factor(Ascites)1                0.441740  1.555411  0.179070  2.467 0.013631
## factor(Ascites)2                0.978996  2.661782  0.181549  5.392 6.95e-08
## Age                             0.048780  1.049989  0.006773  7.202 5.93e-13
## factor(Treatment)1             -0.157295  0.854452  0.377810 -0.416 0.677166
## Prothrombin                    -0.014688  0.985419  0.004080 -3.600 0.000318
## factor(Treatment)1:Prothrombin  0.003450  1.003456  0.005522  0.625 0.532106
##                                   
## factor(Sex)1                   ***
## factor(Ascites)1               *  
## factor(Ascites)2               ***
## Age                            ***
## factor(Treatment)1                
## Prothrombin                    ***
## factor(Treatment)1:Prothrombin    
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                                exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1                      1.6795     0.5954    1.3083    2.1562
## factor(Ascites)1                  1.5554     0.6429    1.0950    2.2094
## factor(Ascites)2                  2.6618     0.3757    1.8648    3.7993
## Age                               1.0500     0.9524    1.0361    1.0640
## factor(Treatment)1                0.8545     1.1703    0.4075    1.7918
## Prothrombin                       0.9854     1.0148    0.9776    0.9933
## factor(Treatment)1:Prothrombin    1.0035     0.9966    0.9927    1.0144
## 
## Concordance= 0.695  (se = 0.017 )
## Likelihood ratio test= 129.9  on 7 df,   p=<2e-16
## Wald test            = 132.4  on 7 df,   p=<2e-16
## Score (logrank) test = 142.5  on 7 df,   p=<2e-16
<<<<<<< HEAD

Opgave 6

=======

Opgave 6

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##                         coef exp(coef)  se(coef)      z        p
## factor(Treatment)1  0.067070  1.069371  0.117868  0.569   0.5693
## factor(Sex)1        0.525289  1.690948  0.126974  4.137 3.52e-05
## factor(Ascites)1    0.441588  1.555174  0.179236  2.464   0.0138
## factor(Ascites)2    0.974063  2.648684  0.181719  5.360 8.31e-08
## Age                 0.049151  1.050378  0.006764  7.266 3.70e-13
## Prothrombin        -0.012923  0.987161  0.002917 -4.431 9.39e-06
## 
## Likelihood ratio test=129.5  on 6 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) + 
##     Age + Prothrombin, data = Cdata)
## 
##   n= 488, number of events= 292 
## 
##                         coef exp(coef)  se(coef)      z Pr(>|z|)    
## factor(Treatment)1  0.067070  1.069371  0.117868  0.569   0.5693    
## factor(Sex)1        0.525289  1.690948  0.126974  4.137 3.52e-05 ***
## factor(Ascites)1    0.441588  1.555174  0.179236  2.464   0.0138 *  
## factor(Ascites)2    0.974063  2.648684  0.181719  5.360 8.31e-08 ***
## Age                 0.049151  1.050378  0.006764  7.266 3.70e-13 ***
## Prothrombin        -0.012923  0.987161  0.002917 -4.431 9.39e-06 ***
## ---
<<<<<<< HEAD
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## 
##                    exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1    1.0694     0.9351    0.8488    1.3473
## factor(Sex)1          1.6909     0.5914    1.3184    2.1688
## factor(Ascites)1      1.5552     0.6430    1.0945    2.2098
## factor(Ascites)2      2.6487     0.3775    1.8550    3.7819
## Age                   1.0504     0.9520    1.0365    1.0644
## Prothrombin           0.9872     1.0130    0.9815    0.9928
## 
## Concordance= 0.696  (se = 0.017 )
## Likelihood ratio test= 129.5  on 6 df,   p=<2e-16
## Wald test            = 132.3  on 6 df,   p=<2e-16
## Score (logrank) test = 142.4  on 6 df,   p=<2e-16
new.age=(65-57)*365

new.surf<-survfit(Cox.fit,newdata = data.frame(Treatment=0,Sex=1,Age=57,Ascites=1,Prothrombin=85))
autoplot(new.surf,xlab="Days",ylab="Overall survival probability")+
  geom_vline(xintercept = new.age)
<<<<<<< HEAD

=======

>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4